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ABSTRACT 


The  Free  Electron  Laser  (FEL)  is  a  candidate  for  a  future  close-in  weapon  system 
that  will  provide  a  longer  protective  range  for  missile  destruction.  The  FEL  is  also 
tunable  to  wavelengths  that  would  give  good  atmospheric  transmission  and  optimal  target 
absorption  characteristics  at  the  target.  This  thesis  describes  single-mode  and  multimode 
simulation  results  of  the  Thomas  Jefferson  National  Accelerator  Facility  (TJNAF)  FEL 
operating  at  far  infrared  wavelengths.  The  TJNAF  FEL  uses  inverse  tapering  and  is 
driven  by  34.5  MeV  and  47.5  MeV  energy  electron  pulses.  Steady-state  power,  weak- 
field  steady  state  gain,  electron  beam  energy  spread  and  optical  spectrum  widths  were 
explored  as  a  function  of  the  desynchronism  and  tapering  rate.  The  simulations  described 
FEL  pulse  evolution  and  short  pulse  effects.  The  simulation  results  have  been  presented 
at  an  International  Conference  held  at  Duke  University,  Durham,  NC  in  August  2000.  In 
addition,  the  results  of  damage  to  Slip-cast  Fused  Silica  samples  by  the  TJNAF  FEL, 
with  and  without  the  effect  of  airflow  are  analyzed.  A  comparison  with  older  damage 
experiments  was  done  in  order  to  develop  scaling  rules  in  the  future. 
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EXECUTIVE  SUMMARY 


Anti-ship  missiles  (ASM)  continue  to  proliferate  and  the  number  of  missiles  that 
come  to  production  using  highly  sophisticated  technology  has  increased  dramatically.  In 
Chapter  I,  through  a  series  of  Phalanx  simulations,  I  show  that  even  though  a  hard  kill  is 
achieved,  the  distance  is  too  short  to  survive  the  missile  debris  that  finally  hits  the  ship 
causing  severe  damage.  Thus,  we  came  to  the  conclusion  that  the  current  close-in 
weapon  systems  (CIWS)  are  inadequate  and  new  means  are  required  to  face  effectively 
the  modem  ASM  threat. 

A  high-energy  laser  CIWS  might  be  a  solution  because  it  provides  an  almost 
instantaneous  reaction  at  the  speed  of  light  and  a  large  missile  destruction  range,  rapid  re¬ 
engagement,  a  large  and  renewable  magazine,  line  of  sight  accuracy,  precision  aim 
pointing  and  single  shot  cost-  effectiveness.  The  Free  Electron  Laser  (FEL)  is  also  a 
candidate  for  CIWS,  because  it  is  tunable  to  various  wavelengths.  The  FEL  can  be  tuned 
to  a  wavelength  that  would  give  good  atmospheric  transmission,  or  is  appropriate  to 
specific  target  absorption  characteristics. 

In  Chapter  II  the  FEL  theory  used  in  this  thesis  is  reviewed  in  detail.  This  theory 
includes  the  pendulum  and  wave  equation  derivations,  the  phase-space  evolution, 
undulator  tapering,  short  pulse  evolution,  and  short  pulse  effects  such  as  the  trapped- 
particle  instability  and  limit-cycle  behavior. 

Chapter  III  presents  new  single  mode  and  multimode  simulation  results  of  the 
TJNAF  FEL  operation  at  far  infrared  wavelengths,  using  inverse  tapering  and  driven  by 
34.5  and  47.5  MeV  energy  electron  pulses.  Steady-state  power,  weak-field  steady-state 
gain,  electron  energy  spread  and  optical  spectrum  width  were  determined  as  a  function  of 


xi 


the  desynchronism  d  and  tapering  rate  S,  using  FEL  pulse  evolution  simulations  and  short 
pulse  effects  described  in  detail  for  each  case.  These  results  are  new. 

In  Chapter  IV,  laser  propagation  through  the  atmosphere  was  reviewed  showing 
that  laser  transmission  is  highly  dependent  on  wavelength.  Then,  we  describe  the  thermal 
diffusion  length  and  showed  that  in  order  to  minimize  thermal  diffusion  and  damage 
sample  materials,  the  laser  spot  size  on  the  target  must  have  diameter  larger  than  the 
thermal  diffusion  length. 

Chapter  V  presents  the  experimental  procedure  and  the  results  of  irradiations 
induced  on  Slip-Cast  Fused  Silica  samples,  with  or  without  the  presence  of  airflow.  A 
comparison  with  older  experiments  was  done  in  order  to  develop  scaling  rules  in  the 
future.  The  first  sample  analyzed  was  irradiated  at  a  wavelength  >1  =  3.10  pm  with  pulse 
repetition  frequency  (PRF)  of  18.7  MHz,  an  average  power  of  105  W,  with  an  average 
flux  intensity  of 490  W/cm2.  The  second  sample  analyzed  had  the  same  irradiation  data 
with  the  first  except  that  the  average  flux  intensity  was  increased  to  10  kW/cm2.  Finally, 
the  third  sample  analyzed  at  the  same  wavelength  with  the  PRF  shifted  to  37.425  MHz 
and  the  average  power  increased  to  500W  for  an  average  flux  intensity  of  10  kW/cm2. 

Chapter  V  includes  the  conclusions  resulting  from  analysis  of  the  new  research 
regarding  the  single  mode  and  multimode  simulations  of  the  TJNAF  FEL  operation,  and 
the  damage  induced  from  the  irradiations  on  the  Slip-cast  Fused  Silica  samples. 

For  FEL  operation  driven  by  the  34.5  MeV  energy  electron  pulses,  it  was  found 
that  even  though  the  maximum  power  is  obtained  for  no  taper  <5  =  0,  negative  taper  of 
3  =  -4n  provides  higher  power  at  desynchronism  values  where  the  operation  is  more 

xii 


stable.  Additionally,  the  maximum  power  with  energy  spread  less  than  6%  is  P  =  165  for 
taper  S  =  -4%  at  desynchronism  d=  0.01 . 

For  FEL  operation  driven  by  the  47.5  MeV  energy  electron  pulses,  it  was  found 
that  negative  tapering  of  <5  =  -4%  gives  the  highest  power  and  stability  in  the  high  power 
operation  area  up  to  d  =  0. 1 6.  The  maximum  power  with  energy  spread  less  than  6%  is 
p  =  180  for  taper  S  =  -4%  at  desynchronism  d=  0.0075.  Keeping  Ay/y  less  than  6%  is  a 
goal  effective  recirculation  of  the  TJNAF  electron  beam. 

In  the  damage  studies,  comparing  our  experimental  results  with  older 
experiments  on  the  same  material,  we  found  that  shifting  the  wavelength  from  X  =  4.825 
pm  to  X  =  3.10  pm  and  the  PRF  from  37.4  MHz  to  18.7  MHz  gave  an  improvement  in 
penetration  rate  of  15.4%  without  airflow  and  25.4%  with  airflow.  In  both  experiments 
the  average  power  was  lOOWatts  and  the  intensity  was  10  kW/cm2.  Shifting  the 
wavelength  from  X  =  4.825  pm  to  X  =  3.10  pm,  increasing  the  power  from  100  Watts  to 
500  Watts,  and  keeping  the  PRF  constant  at  37.4  MHz  with  10  kW/cm2  irradiation 
intensity,  improved  the  penetration  rate  by  1530%  with  the  presence  of  airflow  and  by 
9160%  without  airflow.  The  volume  of  total  damage  increases  by  8.7%  with  the  presence 
of  airflow  and  by  14.1%  without  it,  and  the  volume  of  the  damage  hole  is  increased  by 
525%  with  the  presence  of  airflow  and  by  59%  without  it. 
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I.  INTRODUCTION 

“Zippo  one!  Missiles  inbound  bearing  two  seven  five.”  A  Principal  Warfare 
Officer  or  Tactical  Action  Officer  has  little  time  to  react  to  the  threat  of  an  incoming  anti¬ 
ship  cruise  missile.  He  has  to  decide  within  seconds  about  ship  maneuvers,  and  hard  or 
soft  kill  weapon  deployment.  He  must  act  taking  under  consideration  other  factors  such 
as  threat  type,  wind  conditions  and  his  own  ship  characteristics  (weapon  blind  arcs  and 
mutual  interference).  There  is  really  little  time  for  a  human  operator  to  evaluate  the  threat 
and  respond  accordingly.  Any  error  can  be  fatal.  Anti-ship  missiles  continue  to 
proliferate  and  constitute  a  growing  threat  to  warships.  Advanced  production  technology 
makes  missiles  faster,  stealthier,  lower  flying,  less  susceptible  to  countermeasures  and 
more  agile  in  the  terminal  phase. 

Thus,  we  understand  the  importance  of  a  highly  sophisticated  close-in  weapon 
system  that  provides  quick  reaction  and  a  last-ditch  inner  layer  protection,  against  threats 
that  have  penetrated  through  area  and  point  defenses.  Close-in  weapons  systems  can  be 
autonomous  and  automated.  They  have  their  own  radar  sensors  and  a  multi-barrel  gun 
capable  of  very  high  rates  of  fire,  which  can  tear  apart  the  body  of  an  incoming  missile  at 
ranges  of  about  0.75  to  1.0  nm  (1.4  to  1.8  km).  The  application  of  high-power  lasers  to 
anti-missile  defense  has  been  under  consideration.  Indications  are  that  the  technology 
will  have  matured  to  a  point  where  production  systems  should  be  available  in  the  2010 
time  frame. 
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There  are  six  established  Close-in  weapon  systems  in  service:  Holland’s  Signaal 
Goalkeeper,  Russia’s  Tulamashzavod  AK-630  and  the  Kashtan/Korita  CADS-N-1; 
Spain’s  Meroka,  from  Bazan;  the  Oerlikon-Contraves  Seaguard;  and  the  Raytheon  Mk  15 
better  known  as  Phalanx  which  is  the  most  widely  deployed  CIWS  in  the  world. 

A.  THE  PHALANX  CLOSE-IN  WEAPON  SYSTEM 

1.  Description 

The  Phalanx  CIWS  is  a  closed-loop  weapon  system.  This  means  that  it  is  tracking 
both  the  target  and  the  stream  of  rounds.  It  is  designed  to  provide  the  innermost  layer  of 
defense  against  anti-ship  missiles.  Phalanx  uses  a  Ku-band  pulse  Doppler  sensor  with  a 
tracking  antenna  at  the  front  and  a  search  antenna  at  the  top.  The  system  is  capable  of 
search,  detection,  threat  evaluation,  acquisition,  track  and  fire  modes  which  can  be  done 
automatically  through  a  digital  computer. 

Phalanx  is  normally  set  up  to  scan  a  particular  sector  and  will  automatically 
engage  fast  air  targets  unless  the  hold  fire  button  is  pressed.  Targets  can  be  detected  at 
5.6  km  (3.0  nm)  and  acquired  at  4.3  km  (2.3  nm).  When  the  gun  fires,  the  radar  tracks 
the  centroid  of  outgoing  projectiles,  predicts  their  point  of  closest  approach  to  the  target 
and  corrects  the  aim  of  the  following  projectiles.  This  technique  uses  variable  Pulse 
Repetition  Frequency  (PRF)  with  selected  spectral  frequency  line  tracking  to  measure  the 
stream  of  projectiles’  angular  error.  Firing  usually  begins  at  1 .85  km  (1 .0  nm),  and 
system  reaction  time,  between  target  acquisition  and  fire,  is  reported  to  be  3  seconds. 
Finally  the  gun  caliber  is  20  mm  and  the  rate  of  fire  is  3,000-4,500  rds/min. 
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2.  Phalanx  Simulations 

A  close-in  weapon  system  has  never  been  challenged  in  battle  conditions. 
However  simulations  have  been  used  to  predict  its  performance  against  real  world  threats. 
I  used  these  simulations  [7],  to  explore  the  performance  of  a  typical  close-in  weapon 
system,  the  Phalanx  gun.  Table  1  shows  the  simulations  input  data.  It  must  be  mentioned 
that  the  effects  of  gravity  and  air  drag  have  been  taken  into  account  and  also  that  the 
missile  doesn’t  maneuver. 


Height  of  Gun 

25  m 

Height  of  target 

7m 

Missile  Radius 

0.2  m 

Missile  SDeed 

500  m/s 

Dispersion  angle 

0.002  rad 

Firing  rate 

50  rds/s 

Table  1.  Simulations  Input  Data 

The  Phalanx  gun,  as  in  every  gun  system,  with  a  high  firing  rate,  suffers  from 
dispersion  as  a  result  of  vibrations  shaking  the  barrel.  The  effect  of  the  dispersion, 
although  it  is  small  (1  to  3  millirads),  is  that  the  rounds  do  not  go  precisely  where  they  are 
directed.  Figure  1  shows  the  transverse  positions  of  500  shots  at  a  distance  of  1200  m 
from  the  gun.  As  it  can  be  seen,  the  area  that  that  is  covered  by  the  shots  at  this  distance 
is  about  81  m2  and  only  5  out  of  500  shots  have  hit  the  target  of  radius  0.2m  (shown  as  a 
circle  in  the  center). 
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Figure  1.  Transverse  Positions  of  500  Shots  due  to  Dispersion  at  1200  m 
Range 

Figure  2  shows  the  probability  of  one  shot  hitting  the  missile  vs.  range.  The 
probability  reduces  as  range  increases.  The  probability  is  62.3%  at  a  distance  of  100  m, 
and  drops  to  15.6  %  at  250  m,  and  eventually  becomes  about  0  %  after  1250  m. 
Unfortunately,  one  shot  is  not  enough  to  kill  the  missile.  It  is  found  that  six  to  eight  shots 
are  typically  required  to  either  destroy  the  missile’s  aerodynamics  or  ignite  its  warhead. 
Figure  3  shows  the  number  of  bullets  that  typically  hit  the  missile  as  it  approaches  the 
ship  vs.  range.  It  can  be  seen  that  the  missile  acquired  six  hits  by  the  time  it  reached  a 
distance  of  150  m.  Is  that  far  enough  away  to  say  that  you  have  a  reliable  and  trustworthy 
CIWS  weapon  system? 
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Figure  2.  Single  Round  Intercept  Probability  vs.  the  Range  from  the  Ship 


Figure  3.  Number  of  Accumulated  Hits  vs.  the  Range  from  the  Ship 


The  answer  is  “No.”  Because  even  if  you  kill  the  missile  at  that  distance  it 
doesn’t  mean  that  you  survived  the  attack.  The  debris  fuel,  and  warhead,  from  the 
destroyed  missile  can  still  hit  the  ship  causing  severe  damage.  Figure  4  shows  the 
random  trajectories  of  20  random  sized  debris  fragments  coming  towards  the  ship,  from  a 
missile  explosion  at  190  m  from  the  ship.  As  it  can  be  seen,  10  out  of  20  sample 
fragments  have  hit  the  ship.  Figure  5  shows  the  probability  of  a  fragment  hitting  the  ship 
vs.  range.  At  the  missile  killing  distance  of  150  m  found  earlier  60  %  of  the  missile 
debris  hit  the  ship. 


Figure  4.  Possible  Debris  Trajectories  from  a  Destroyed  Inbound  Missile 
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Figure  5.  Probability  of  a  Fragment  Hitting  the  Ship  vs.  Destruction  Range 

The  scenario  explored  at  the  previous  simulations  was  simple.  The  missile  was 
non-maneuvering,  the  engagement  commenced  at  maximum  range,  the  threat  was  not 
multiple  (one  missile  was  encountered)  with  perfect  tracking  by  the  CIWS  radar. 
However,  the  results  are  very  disappointing  and  the  existing  CIWS  has  been  proven 
inadequate.  A  new  CIWS  system  is  required,  and  therefore  research  should  be  done  on  a 
point  air  defense  system  that  is  much  more  accurate  with  longer  range. 

B.  HIGH  ENERGY  LASERS 

Lasers  are  distinguished  by  the  lasing  medium  they  use.  Chemical  lasers  have 
become  the  laser  for  military  directed  energy  applications.  The  lasing  mediums  used  are 
usually  hydrogen  fluoride  (HF),  deuterium  fluoride  (DF),  and  oxygen  iodine.  For 
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example,  the  MIRACL  (Mid  Infrared  Advanced  Chemical  Laser),  is  a  2.2  MW  deuterium 
fluoride  laser,  with  laser  wavelength  3.800  nm,  and  has  successfully  downed  several 
dummy  missiles  at  a  half  mile  range  in  test  demonstrations  at  the  White  Sands  Missile 
Test  Range  in  the  USA. 

The  laser  device  transforms  energy  from  a  power  source  into  light  energy  of  a 
single  wavelength  propagating  as  a  highly  directional  beam.  Every  laser  oscillator  has 
four  components:  an  amplifying  medium  that  increases  the  intensity  of  a  light  beam  as  it 
passes  through  the  medium,  an  energy  source  that  pumps  the  electrons  in  the  amplifying 
medium  to  higher  energy  levels,  an  optical  cavity  consisting  of  an  arrangement  of  mirrors 
that  feed  the  amplifying  medium,  and  finally  a  mechanism  that  leaks  a  useful  part  of  the 
beam  of  light  from  the  optical  cavity  without  significantly  diminishing  the  amplification 
process  that  takes  place  in  the  cavity  and  laser  medium.  The  effect  of  this  mechanism  is 
to  provide  a  beam  of  coherent  light  energy,  which  means  that  the  emitted  light  energy  is 
of  the  same  frequency  and  is  in  phase. 

A  laser  “kills”  by  heating  a  target’s  surface,  possibly  at  a  vulnerable  point  such  as 
a  missile’s  solid  rocket  casing  or  liquid  fuel  tank.  The  laser  must  fulfill  several  tasks  to 
kill  a  target.  First  it  must  generate  a  high-energy  laser  beam.  Then,  the  acquisition, 
tracking  and  pointing  system  acquires  and  tracks  and  hands  over  the  aim  point  to  the 
beam  control  system  and  range  information.  Finally  the  beam  control  system  focuses  the 
beam  on  the  target,  and  maintains  the  laser  spot  on  the  same  area  of  the  target. 

Developing  a  HEL  will  give  to  the  ship’s  defense  a  potent  capability.  HEL 
weapons  will  be  able  to  deliver  destructive  amounts  of  energy  to  a  target  up  to  10 
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kilometers  away  at  the  speed  of  light,  with  enough  energy  to  destroy  a  ballistic  missile  in 
a  few  seconds.  The  accuracy  needed  is  tens  of  nanoradians  in  angle.  It  will  be  able  to  kill 
a  target  in  a  few  seconds  and  shift  to  a  new  one.  Unlike  guided  missiles,  the  HEL’s 
sophisticated  and  expensive  computing  and  mechanical  systems,  which  hold  laser  energy 
accurately  onto  the  target  are  on  the  launch  platform  and  are  not  destroyed  during 
engagement. 

Concluding,  the  laser  offers  many  advantages  over  a  gun  and  missile-based  close- 
in  weapon  systems.  These  include  almost  instantaneous  reaction  at  the  speed  of  light  and 
long  missile  destruction  ranges,  rapid  re-engagement,  a  large  and  renewable  magazine, 
line  of  sight  accuracy,  precision  aim  pointing,  and  high  single  shot  cost-  effectiveness. 

In  practice,  various  problems  are  yet  to  be  solved,  such  as  thermal  blooming, 
atmospheric  attenuation,  beam  focusing,  power  generation,  and  equipment  weight  and 
volume.  However,  while  the  technology  promises  much,  research  and  development 
efforts  to  produce  a  trustworthy  and  deployable  system  are  continuing 

C.  THE  FREE  ELECTRON  LASER 

The  free  electron  laser  is  a  high-energy  laser  and  consists  of  three  main 
components:  the  electron  accelerator,  the  undulator,  and  the  optical  resonator.  An 
electron  beam  is  accelerated  in  the  accelerator  to  a  relativistic  speed  and  enters  into  the 
resonator.  Inside  the  resonator  there  is  the  undulator.  The  undulator  consists  of  a  series 
of  alternating  magnets  that  produce  a  transverse,  sinusoidal  magnetic  field  along  the 
undulator  axis,  which  makes  the  relativistic  electron  beam  wiggle.  In  this  way,  the 
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electrons  accelerate  and  radiate  energy  in  the  forward  direction.  In  the  presence  of  light, 
some  electrons  gain  kinetic  energy  from  the  optical  field  and  move  faster,  while  other 
electrons  are  loosing  energy,  giving  it  to  the  field,  and  slowing  down.  This  results  in 
bunching  of  the  electrons,  leading  to  coherent  radiation.  Bunching  is  essential  because 
the  emission  rate  for  the  bunched  electron  beam  is  proportional  to  the  square  of  the 
number  of  electrons.  On  the  other  hand,  the  emission  rate  for  the  beam  of  randomly 
positioned  electrons  is  only  proportional  to  the  number  of  electrons.  Each  electron 
typically  radiates  millions  of  coherent  photons  in  one  pass  through  the  undulator.  Figure 
6  shows  a  typical  configuration  of  an  FEL  oscillator. 


Figure  6.  Oscillating  FEL  Configuration,  From  Ref.  [3] 

The  FEL  has  potential  advantages  over  the  conventional  CIWS  that  are  in  use 
now.  The  first  is  longer  destruction  range,  which  is  essential  for  the  ship’s  survivability 
because  it  prevents  ship  damage  from  debris.  Also,  it  improves  our  ability,  by  giving  us 
time,  to  engage  a  second  and  third  successive  incoming  missile.  The  second  is  the 
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precision  in  pointing  and  shot  accuracy  because  lasers  don’t  suffer  from  dispersion,  and 
can  focus  on  an  incoming  missile  at  long  ranges. 

An  important  advantage  the  FEL  has  over  the  other  laser  weapons  is  that  it  is 
tunable  to  specific  wavelengths.  The  irradiation  wavelength  of  the  FEL  can  be  tuned  in 
three  ways,  either  by  changing  the  undulator  wavelength,  or  by  changing  the  undulator’s 
magnetic  field  amplitude,  or  by  changing  the  accelerator’s  electron  beam  energy.  The 
tunability  gives  advantage  in  selecting  a  wavelength  that  would  give  good  atmospheric 
transmission  and  is  appropriate  to  specific  target  absorption  characteristics. 

D.  PURPOSE  OF  THIS  THESIS 

This  thesis  consists  mainly  of  two  parts.  In  the  first  part,  multimode  simulations 
are  used  to  analyze  the  operation  of  the  Tomas  Jefferson  National  Accelerator  Facility 
(TJNAF)  FEL  with  a  negatively  tapered  undulator.  The  evolution  of  short  optical  pulses 
in  the  far  infrared  are  described  and  show  the  effects  of  taper  on  gain,  power  and 
desynchronism.  Additionally,  single-mode  simulations  are  used  to  explore  the  effects  of 
negative  taper  on  gain.  In  the  second  part,  the  results  of  the  TJNAF  FEL  irradiations  on 
samples  of  Slip-cast  fused  Silica  are  described  and  discussed. 
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II.  THE  FREE  ELECTRON  LASER  THEORY 


A.  THE  PENDULUM  EQUATION 

When  the  electron  beam  enters  in  the  undulator,  the  interaction  between  the 

electron  beam  and  the  optical  field  takes  place  along  the  undulator.  A  Cartesian 
coordinate  system  is  set  at  the  beginning  of  the  undulator  with  the  z  direction  along  the 
undulator  axis.  Assume  a  helical  undulator  magnetic  field  B : 

B  =  (Bx,By,Bz)  =  B(cos(k0z),sm(k0z),0)  (1) 

where  k0  =  2tt/  is  the  undulator  wavenumber,  z  is  the  distance  down  the  undulator, 
and  B  is  the  magnetic  field  strength  and  Bx,  By,  and  B,  are  the  magnetic  field 
components  in  the  x,  y  and  z  direction.  Assume  a  corresponding  circularly-polarized 
plane  wave  for  the  optical  field  that  is  generated  inside  the  undulator 

Es  =  E{  cos('F),-sin(T'),0),  (2) 

Bs  =E(sin('P),cos(vP),0),  (3) 

where  E  is  the  electric  and  magnetic  field  amplitude  in  cgs  units,  Y  =  kz-cot  +  <f>,  <j>  is 
the  optical  phase  angle,  co  is  the  angular  frequency,  t  is  the  time  and  k  =  In  /  X  is  the 
optical  wave  number.  The  relativistic  Lorentz  force  equations  that  govern  the  motion  of 
the  electron  in  the  undulator  are 

=  -—[Es  +  fi  x  (Bs  +  B)l  (4) 

dt  me 
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(5) 


dy 

~dt 


me 


y~2=\-pp,  (6) 

where  m  is  the  electron  rest  mass,  e  is  the  electron  charge  magnitude,  ft  =  v!  c  is  the 
dimensionless  electron  velocity,  and  y  =  (1  -  p  ■  pyvl  is  the  Lorentz  factor. 

Substituting  the  fields  in  Equation  (1),  (2)  and  (3)  into  Equation  (4)  and  (5),  and 
split  p  into  two  components:  PL=(px,py,  0)  and  P2  =(0,0  ,pz)  so  that 

$  ( yR  }  £ 

= - W  -  Pz  )(cos(T),-sin(vF),0)  +  BPz(-sm(k0z),cos(k0z),0)],  (7) 


d  ( yB  )  e 

dT  =  ~~^[E{Px  C0S(VF)  ■  p> sin(T))  +  B{f3*  siri(k°z)  -  P y  cos(^oz))]>  (8) 

and  ~  =  ~—E[px  COS(T)  -  P  sin('E)].  (9) 

at  me 

For  relativistic  electron  (y  » 1),  E(  1  -  p.)  «  Bp.  and  Equation  (7)  becomes 

dfyB, )  e 

— T  = - 5^z[-sm(A:0z),cos(A:0z),0].  (10) 

dt  me 

Equation  (10)  can  now  be  integrated  with  respect  to  t  to  find  that 

-  K 

Px= - [cos(£0z),sin(£0z),0].  (11) 

r 

where  the  undulator  parameter  is  K  =  eBXa  Hume1.  Inserting  Equation  (1 1)  into 
Equation  (9)  we  find  that 
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(12) 


dy  eEK 
dt  ymc 

where  the  electron  phase,  C,  -  (k  +  k0)z  -  cot, 


■cos  (£  +  0), 

determines  the  z  position  relative  to  an 


optical  wavelength,  X 

Substituting  Equation  (1 1)  in  Equation  (6)  we  obtain 

y-2  =1-P-P  =  \-P2z-Pl=\-Pl  ~K2/y2,  or  (1  +  K2)y~2  =1  -p2.  (13) 


Differentiating  both  sides  of  Equation  (13)  with  respect  to  t ,  dy/dt  and  d(3-/dt  are  related 


f  _  r2PzPz 
7  d  +  ^2)' 


(14) 


Now  taking  the  first  and  the  second  time  derivative  of  the  electron  phase,  £  we  obtain 


£  =  (k  +  k0)cp,  -a. 

(15) 

C=(k  +  k0)cpz. 

(16) 

Solving  Equation  (16)  for  fi2  and  substituting  in  Equation  (14),  we  get 

r 

y  (1  +  K2)(k  +  k0)c' 

For  relativistic  electrons  pz  =  1  and  k»k0,  Equation  (17)  becomes 

t=A L 

7  2u>0 

where  near  resonance  co0  =  a>(\  +  K 2 )  /  2y 2 .  The  resonance  condition  is 


(17) 


(18) 


Z  = 


a0(i  +  a:2) 
2/ 2 


(19) 


In  order  for  the  electrons  to  interact  strongly  with  the  optical  field,  they  must  satisfy  the 
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resonance  condition:  namely  that  one  wavelength  of  light,  \  passes  over  an  electron  as 
the  electron  travels  through  one  undulator  wavelength,  V 


Solving  Equation  (18)  for  C,  and  substituting  into  Equation  (12)  we  get  the 
Pendulum  Equation  of  electron  motion  within  the  undulator 

2  coj  leKEco, 

$  = - = - 5 - cos(£  +  $). 


r 


y2mc 


(20) 


It  is  helpful  when  designing  a  FEL  to  use  dimensionless  parameters  when  ever  possible. 
Therefore,  we  define  the  dimensionless  time,  r  =  ct  /  L,  where  L  is  the  undulator  length, 
so  that  ris  equal  to  0  at  the  beginning  of  the  undulator  and  1  at  the  end.  We  also  define 
the  dimensionless  complex  optical  field  a  =  \a\e‘*  where  |a|  =  4neNKL\E\/ y2mc2 .  Thus 
the  Pendulum  Equation  becomes 


v  ~  C  =|a|c°s(4'  +  ^).  (21) 


o 

where  (...)=  d{.)/  dr,  and  the  dimensionless  electron  phase  velocity  is 


v  —  <->  -  L[(k  +  K )Pz  —  &];  if  v—0  the  electrons  are  at  resonance.  The  electron  phase 
velocity  is  proportional  to  the  electron  dimensionless  velocity  j3z,  so  that  when  v  decreases 
the  electron  slows  down  and  gives  energy  up  to  the  optical  field.  This  is  what  we  always 
want  to  occur  for  a  powerful  laser.  In  order  to  have  optical  gain,  more  electrons  must 
give  up  energy  to  the  optical  field  than  take  energy  from  it.  From  the  Pendulum 
Equation,  we  see  that  the  maximum  energy  loss  from  the  electrons  occurs  when 
cos(^  +  0)  =  -1  or  when  electrons  are  bunched  at  £  +  <f>  =  n.  This  is  the  goal  of  the  FEL. 
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B.  THE  WAVE  EQUATION 

The  development  of  the  complex  optical  wave  equation  starts  from  Maxwell’s 
wave  equation  acting  on  the  vector  potential  A , 


_Lil 

c2  dt 2 


-*  4/r  - 

A  =  -—Jx, 
c 


(22) 


where  J±  is  the  transverse  current  density  and  the  vector  potential  corresponding  to  the 
fields  in  Equation  (2)  and  (3)  is 

A(z,t)  =  -£(z,t)[sm('¥),cos('¥),0].  (23) 

co 

where  optical  phase  is  'F  =  kz  -  cot  +  <f)(z,t).  Taking  the  second  derivative  of  Equation 
(23)  with  respect  to  z  gives 

zr  =  ttt[* + itIM'1')  -siitfno] 

dz  k  8z  L  dz  _ 

+-^[sinCP),  cos(T'),0]+ ~\k  +  |^l[cos(T)-sin(^),0] 
k  dz  k  dz\_  dz 

+  — ^[cos(T')-sin(T'),0]+  — \k  +  ^-]  [-sin(T)  -cosCno]  (24) 
k  dz  k  L  dz  J 

Now  taking  the  second  derivative  of  Equation  (23)  with  respect  to  t  gives 

— =  ^T\^~~co  [cos(lF) -sin(lF),0]+-^-f-[sin(xF),cos(T,),0] 
c  dr  dt  L  dt  J  dt 

^  -  6>l[cos(lF)-sin(xf/),0]+ E  Jcos('F),-  sin('F),0] 
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+  E  ^--co  [-sin^-cos^XO]  (25) 

ot  J 

Next,  assume  the  optical  amplitudes  and  phases  varying  slowly  in  time  and  space 
so  that  dE  /  dz  «  kE,  d(f)  /  dz  «  k(j>,  8E  /  dt «  coE,  d<f>!  dt  «oxf>,  and  w  =kc.  This 
allows  us  to  use  only  the  first-order  derivatives  because  the  second  order  derivatives  are 
small.  Applying  this  approximation  to  Equations  (24)  and  (25)  and  inserting  the  results 
into  Equation  (22)  we  obtain 

1  32  l-  JdE  \dd)\  -j 

{&+cars(nsin(no] 

+  2£  [- sin('F),-cos('i'),0]  (26) 

The  current  density  for  a  single  electron  is  J±  =  -ecJ3±.  Substituting  Equation 
(11)  for  /?x,  we  get 


J±= - [cos(£0z),  sin(&0z),0}  (27) 

7 

We  introduce  the  dimensionless  time  r  =  ctIL  again  so  that  r  is  equal  to  0  at  the 
beginning  of  the  undulator  and  1  at  the  end.  Substituting  r  and  Equation  (27)  in 
Equation  (26),  we  obtain 

2  TIT  [cos(^)>sin(y),o]+2£  [-sin('F),-cos(vf'),0]s 

_L  OT  J  _Z,  OT  J 
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4  mK 


[cos(£0z),sin(£0z},0}. 


Equation  (28)  can  be  separated  using  two  expressions  representing  change  in  E  and 
change  in  0 .  Thus  Equation  (28)  becomes 


dE  IneLK 


cos{£  +  0), 


50  _  2neLK 
5t  y 


sin(£  +  0), 


where  £  =  (k  +  k0)z  -  cot  is  the  electron  phase.  Equations  (29)  and  (30)  are  the  wave 

equations  driven  by  a  single  electron,  but  we  need  to  evaluate  the  result  from  a  beam  of 
electrons  using  a  sum  over  many  electrons  in  the  FEL  beam.  This  can  be  done  by  taking 
the  average  <...>  and  multiplying  it  by  electron  density,  p ,  resulting  in 


dE  InepLK 


<  cos^  +  0)  >, 


50  2mpLK 

~Z~  = - - - <  sm(<^  +  0)  >  ■ 

or  y 


Taking  the  derivative  of  the  complex  electric  field  E  =  \E\e'*  with  respect  to  t 


and  inserting  Equations  (31)  and  (32)  we  obtain 


—  =  [<  cos(^  +  0)  >  —i  <  sin((^  +  0)  , 

5t  Y 


which  can  be  simplified  to 
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(34) 


dE_  =  _lnepLK_  _K 


dr 


<  e  *  >. 


From  the  pendulum  equation  we  know  that  the  dimensionless  optical  field  is  a  =  Ja|e(> , 
where  |a|  =  AweNKL\E\l  y~mc~ .  Therefore  the  final  form  of  the  wave  equation  is 


da  ° 

—  =  a  =  -j  <  e~'^  >, 
dr 


(35) 


where  j  =  8;r  e~  pNK2  L2  / y3mc2  is  the  dimensionless  current. 

The  growth  of  the  optical  field  and  therefore  the  development  of  the  laser  energy 
depend  on  both  the  dimensionless  current  and  the  average  electron  phase.  If  there  is  no 
current  or  no  electron  bunching,  there  is  no  growth  of  the  optical  field. 


C.  THE  FEL  PHASE-SPACE  EVOLUTION 

The  pendulum  equation  motivates  a  discussion  of  the  FEL  electron  phase-space 
evolution.  The  phase-space  evolution  describes  the  microscopic  motion  of  the  electrons 
on  the  scale  of  the  optical  wavelength.  The  phase-space  simulations  that  will  be  shown 
are  at  the  upper  limit  of  the  weak-field  regime  (a0  =  3)  and  low  current  (j  =  1). 

Figure  7  shows  the  FEL  phase-space  evolution  of  twenty  sample  electrons  at 
resonance.  There  are  typically  many  millions  of  electrons  in  each  optical  wavelength,  but 
the  evolution  of  only  twenty  sample  electrons  are  shown.  The  vertical  axis  in  phase- 
space  is  proportional  to  the  electron  energy  while  the  horizontal  axis  is  the  electron’s 
phase  within  a  section  of  the  electron  beam  one  wavelength  of  light  long.  Figure  7  also 
shows  the  “separatix.”  The  separatix  is  given  by  the  formula  v  =  2|a|[l  +  sin(^  +  <p)]  [8], 
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It  separates  the  open  and  closed  phase-space  paths  and  passes  through  the  critical  points 

I  |l/2 

(-x/2,0)  and  (3tr/2,0).  The  peak-to-peak  height  of  the  separatix  is  4|a|  ‘[8].  On  the  top 

right  of  the  phase-space  evolution  the  gain  in  optical  energy  is  plotted,  and  on  the  right 
bottom  the  change  of  the  optical  field  phase.  Starting  at  resonance  leads  to  no  gain.  This 
happens  because  the  half  of  the  electrons  gain  energy  from  the  optical  field,  while  the 
other  half  give  away  the  same  amount  of  energy  to  it.  It  can  also  be  seen  that  at  resonance 
the  optical  phase  shift  is  relatively  large  and  the  electrons  become  bunched  at  £  =  tr/2. 


-%/2  C  3 tr/2  0  x  1 

Figure  7.  The  FEL  Phase-Space  Evolution  at  Resonance 

Figure  8  shows  the  gain  and  phase  curves  with  respect  to  the  initial  electron  phase 
velocity  Vq .  The  gain  spectrum  is  anti-symmetric  in  Vq  with  a  peak  gain  of  G  =  0.125  at 
v0  =  2.7.  At  resonance,  Vo=  0,  there  is  no  gain  at  any  time,  while  for  values  0<Vo<6.2 
there  is  net  absorption  of  the  optical  power.  The  optical  phase  shift  is  symmetric  in  Vq 
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with  a  peak  value  of  A<f>  =  0.074  at  Vq  =  0.  At  the  phase  velocity  for  maximum  gain 
Vo  =  2.7,  the  phase  shift  is  only  A (f>  =  0.015 .  Figure  9  shows  the  FEL  phase-space 
evolution  slightly  off  resonance  at  the  optimum  initial  phase  velocity  of  v0  =  2.7.  It  can  be 
seen  that  the  electron  bunching  occurs  at  f =  j,  the  phase  for  which  the  maximum 
electron’s  energy  loss  occurs. 
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Figure  8.  The  FEL  Gain  and  Phase  Curves 
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Figure  9.  The  FEL  Phase-Space  Evolution  at  Optimum  Initial  Phase 
Velocity 


D.  THE  TAPERED  UNDULATOR 

The  FEL  resonance  condition  found  during  the  derivation  of  the  pendulum 
equation  is 

A  =  M±£!)  (36) 

2r2 

where  Xis  the  optical  wavelength,  X?  is  the  undulator  wavelength,  7  is  the  Lorentz  factor 
and  K  =  eBl0 1  lime1  is  the  dimensionless  undulator  parameter.  At  normal  saturation, 

the  electron  beam  loses  enough  energy  to  the  laser  light  to  shift  across  the  gain  bandwidth 
by  A  v  =  AnNAy  /  y  [9],  moving  away  from  the  value  of  v  for  peak  gain  to  a  value  of  v 

for  no  gain  or  even  loss. 
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Taper  the  undulator  can  be  used  to  accelerate  the  electron’s  phase  v  =  C,  =  8  so 
that  resonance  is  maintained.  This  can  be  done  either  by  increasing  (negative  tapering)  or 
decreasing  (positive  tapering)  the  undulator  wavelength,  or  the  undulator  magnetic  field 
strength,  or  both.  The  phase  acceleration  6  is  given  by  8  =  -2;zNAA0  /  A0  when  the 

undulator  wavelength  is  modified,  or  by  8  «  -AtzNK 2  AK  /[(l  +  K 2  )K0  ]  when  the 

undulator  field  strength  is  modified  [9].  Now  that  the  tapering  is  introduced,  the 
pendulum  equation  includes  an  additional  constant  torque  8  and  has  the  form 

oo 

^  =  £  +  |a|cos(^  +  ^).  (37) 

The  desirable  criteria  for  positive  tapering  are  [9] : 

|a|  >  <i>  >  4|a|'/2  >  2rc.  (38) 

The  condition  on  the  left  is  the  requirement  to  have  trapped  electrons  in  the  tapered  phase 
space;  the  middle  requirement  dictates  that  the  tapered  acceleration  exceed  the  untapered 
deceleration;  and  finally,  the  right  condition  dictates  that  the  FEL  work  in  the  strong  field 
regime  (a  »  7r). 

In  the  next  chapter,  single-mode  and  short  pulse  multi-mode  simulations  on 
negative  tapering,  (increasing  the  undulator  magnetic  field),  will  be  discussed.  It  is 
proven  here  that  a  constant  torque  8  corresponds  to  an  approximately  linear  increase  of 
the  magnetic  field  along  the  undulator  axis.  The  electron  Equation  with  tapering  in  weak 
field  optical  fields  has  the  form 
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v  =  S.  (3 

Integrating  Equation  (39)  with  respect  to  dimensionless  time  t,  and  substituting  r=ct/L 
and  v  =  L[(k  +  k0)fiz-k\  we  get 

i[(*+*„)A  -*]=*.  << 

where  I  is  the  undulator  length,  k  is  the  optical  wavenumber,  k0  is  the  undulator 
wavenumber,  and  z  is  the  distance  down  the  undulator.  Inserting  =  1  -  (1  +  K2)/  2y 2 


in  Equation  (40)  we  obtain 


r  ,(,  l  +  K2)  ,  f,  l  +  K2  )  ,  s 

I  2r2  J  °{  2 r2  ). 


Assuming  that  k»  k0  for  y  «  70  and  solving  for  AT,  equation  (41)  becomes 


K( r)  =  . 


*.-^1-1. 


The  resonance  condition  for  the  wavenumbers  is  k  =  k02y2  /(I  +  Kq  ),  where  K0  is  the 
undulator  parameter  in  the  beginning  of  the  undulator.  Solving  the  resonance  condition 
for  2y/k,  and  substituting  to  Equation  (42)  we  get 


Recall  that  k0  =  2n  /  Z0  and  L  =  Ns l0 ,  where  N  is  the  number  of  undulator  periods,  so 


that  Equation  (43)  becomes 
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rl±Kl  ^ 

V2  *NK0J 


5t. 


(44) 


Figure  10  shows  Equation’s  (44)  graph  for  K0=  1,  N=  41,  and  8  =  -8x  For  N  large,  the 
square-root  argument  can  be  accurately  expanded  to 


K(t)  *  K0 


i  .hiU 


\ 


2nNK, 


(45) 


o  J 


which  shows  that  a  straight  line  is  approximately  correct. 
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Figure  10.  Linear  FEL  Inverse  Tapering 

Figure  1 1  shows  the  final  FEL  phase-space  distribution  for  a  periodic  undulator 
(5  =  0),  while  Figures  12  and  13  show  the  final  phase-space  distribution  of  a  positive 
(5  =  47 r)  and  a  negative  (8  =  -4x)  tapered  undulator,  respectively. 
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Figure  13.  The  FEL  Phase-Space  Evolution  of  a  Negative  Tapered 
Undulator  (6  =  -4ir)  in  Strong  Field 

The  electrons  in  the  simulations  started  at  the  optimum  phase  velocity  for  each 
case  and  with  a  Gaussian  spread  of  in  phase  velocity  of  oQ  =1.3.  The  initial  strong 

optical  field  a0  =  40  is  amplified  by  a  current  of  j  =  7,  in  an  undulator  with  N=  41 

periods.  The  separatix  in  the  tapered  undulators,  Figures  12  and  13,  surrounds  a  smaller 
area  of  closed  orbits  than  does  the  untapered  case  in  Figure  1 1 ,  and  is  given  by 

v2  =2S(C -C0)  +  2\a\(sin(c  +  ^)-sin(4'0 +^))  (46) 

where  4"0  =  2n  -  cos  1  (-  5  /|a|)- <f>  [9].  It  must  also  be  noted  that  the  separatix  shapes  of 

the  positive  and  negative  taper  are  reversed;  the  area  enclosed  decreases  as  the  taper 
phase  acceleration  5  increases,  and  increases  in  stronger  fields  as  a0  increases 
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In  the  untapered  undulator  of  Figure  1 1,  the  electrons  in  the  presence  of  strong 
optical  field  have  over  bunched  and  the  gain  is  saturated.  The  gain  at  the  end  of  the 
evolution  is  Gf  =  3.74%.  Using  a  positively  tapered  undulator  (5  =  47t)  of  Figure  12, 

there  is  bunching  of  the  trapped  electrons  in  closed  orbits  inside  the  separatix  near  f  =  ir, 
and  the  final  gain  is  increased  to  Gf  =  5.36%.  For  the  negatively  tapered  undulator 

(5  =  -4tt)  in  Figure  13,  no  electrons  are  trapped  inside  the  seperatix  and  all  of  them  are  in 
open  orbits.  The  bunching  occurs  as  electrons  travel  around  the  closed  orbit  region  of 
phase  space  at  final  phase-space  sites  around  v  «  -12  where  the  electrons  have  lost 
significant  energy,  and  the  final  gain  is  Gf  =  8.45%. 

Concluding  it  must  be  mentioned  that  the  tapered  undulators,  both  positive  or 
negative,  are  more  efficient  in  strong  optical  fields,  but  have  smaller  gain  than  the 
untapered  undulator  in  weak  fields. 

E.  THE  SHORT  PULSE  EVOLUTION 

The  majority  of  the  FEL  oscillators  are  driven  by  short  electron  pulses  rather  than 

a  continuous  beam.  Electron  pulses  are  considered  to  be  short  when  they  are  comparable 
to  one  slippage  distance,  N\  where  N  is  the  number  of  undulator  periods  and  Xis  the 
optical  wavelength. 

As  short  electron  pulses  enter  the  FEL  oscillator,  short  optical  pulses  are 
generated  due  to  spontaneous  emission  bouncing  back  and  forth  between  the  resonator 
mirrors.  The  distance  between  the  mirrors  is  S,  and  is  bigger  than  the  undulator  length, 

S  >  L.  The  optical  pulse  travels  a  distance  25,  and  at  time  intervals  of  25  /  c  reaching 
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the  beginning  of  the  undulator  at  r  =  0.  At  this  time,  an  electron  pulse  also  enters  the 
undulator.  The  desynchronism,  d,  is  as  the  displacement  between  the  electron  pulse  and 
the  optical  pulse  at  r  =  0,  normalized  to  the  slippage  distance.  If  d  -  0,  then  the  optical 
pulse  coincides  with  the  electron  pulse  each  pass;  this  is  called  exact  syncronism. 
Desynchronism  is  easily  adjusted  in  practice  by  moving  the  resonator  mirrors  by  microns 
with  a  piezoelectric  crystal. 

Figure  14  shows  a  pulse  evolution  simulation  with  electrons  entering  in  the 
undulator  at  the  optimum  phase  velocity  v0,  at  exact  synchronism  (d  =  0).  The  peak 

current  is  j  =  6,  the  pulse  length  is  a,  - 1 .3,  the  desynchronism  value  is  d  -  0,  the 
quality  factor  is  0  =  10  implying  1 0%  resonator  loss  per  pass,  the  tapering  phase 
acceleration  is  5  =  -6tt,  the  number  of  undulator  periods  is  N  =  41,  and  finally  the 
standard  deviation  of  a  small  fluctuation  in  the  initial  electron  phases  is  8C,  =  10"4.  The 
simulation  plot  contains  six  windows.  At  the  upper  left  is  the  dimensionless  optical  field 
shape  a (z,n),  as  a  function  of  z  which  is  scaled  to  the  slippage  distance,  and  its 
evolution  with  the  number  of  passes  n.  At  the  upper  middle  and  right  are  the  optical 
power  spectrum  P(v,n)  and  the  electrons’  spectrum  f(v,n)  respectively,  as  a  function  of  the 
phase  velocity  v  and  their  evolution  with  n.  At  the  lower  left,  the  longitudinal  profile  of 
the  current  density  j  is  shown  for  reference  at  times  r  =  0  and  r  =  1,  as  a  function  of  z 
which  is  also  scaled  to  the  slippage  distance.  The  electron  pulse  is  assumed  to  be 
parabolic  in  shape  with  the  form  j(z)  =  y(l  -  2z2 1  a\ )  for  j(z)  >  0,  and  0  otherwise  [10]. 
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At  the  lower  middle  is  the  weak  field  gain  G(v )  for  reference,  as  a  function  of  v.  At  the 
lower  right  is  the  evolution  of  the  total  power  P(n)  as  the  number  of  passes  n  increases. 
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Figure  14.  The  FEL  Pulse  Evolution  at  Exact  Synchronism  (</=0) 

Surprisingly  we  see  in  Figure  14  that  exact  synchronism  leads  to  no  steady  state 
power.  This  result  can  be  explained  in  the  following  way.  The  electron  pulse  and  the 
optical  pulse  enter  in  the  undulator  at  the  same  moment  each  pass.  They  begin  to  interact 
along  the  undulator,  but  in  the  beginning  no  electron  bunching  occurs,  and  therefore, 
there  is  no  gain.  This  delay  in  the  gain  is  called  lethargy  [10].  Then  bunching  occurs  and 
gain  developes.  As  the  pulses  travel  together,  the  gain  increases  but  since  the  light  pulse 
moves  slightly  faster  than  the  electron  pulse,  the  electron  pulse  falls  behind  most  of  the 
amplification  occurs  at  the  trailing  edge  of  the  light  pulse.  Consequently,  the  light  pulse 
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is  distorted  on  each  pass  and  its  centroid  is  traveling  with  speed  lower  than  c  even  though 
it  is  in  vacuum.  Eventually,  over  many  passes  the  coupling  between  the  two  pulses  will 
be  lost  and  the  resonator  losses  will  cause  the  optical  pulse  to  decay  and  the  FEL  will  fail 
to  operate.  The  simulation  in  Figure  14  has  been  artificially  started  with  a  coherent  pulse 
of  light  and  shows  that  peak  power  occurs  at  n  =  50  passes,  but  then  the  optical  and 
electron  pulses  decouple  so  that  the  FEL  power  dies  out  after  about  n  - 170  passes. 

In  order  to  compensate  for  the  slower  speed  of  the  distorting  light  pulse,  the 
resonator  path  S  must  be  reduced  by  AS  so  that  d  -  -2  AS  /  NX.  If  d  becomes  too  big, 
the  compensation  is  too  severe,  and  the  pulses  decouple  after  many  passes  leading  again 
to  no  steady  state  power.  Figure  15  shows  the  pulse  evolution  simulation  for  an  FEL  with 
exactly  the  same  characteristics  as  that  of  Figure  14,  but  with  desynchronism  d  =  0.015. 
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Figure  15.  The  FEL  Pulse  Evolution  at  <f=0.015 
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Figure  15  also  shows  that  no  power  develops  until  n  =  250  passes,  and  that  steady  state 
power  is  reached  after  n  =  500  passes. 


F.  THE  TRAPPED  PARTICLE  INSTABILITY 

The  trapped-particle  instability  is  an  effect  that  occurs  in  strong  optical  fields 

when  electrons  become  trapped  in  deep  potential  wells  in  phase-space.  The  height  of  the 
separatix  and  area  enclosed  increases  with  the  optical  field  amplitude  a.  When  the  FEL 
reaches  high-power,  the  optical  field  gets  longer  so  that  many  electrons  become  trapped 
on  the  closed  orbits  of  phase-space.  The  electrons  that  are  trapped  near  the  stable  fixed 
point,  electron  phase  C,  «  n  /  2  -  <j>,  are  executing  a  part  of  a  synchrotron  oscillation  [11]. 

The  stable  fixed  point  for  the  tapered  undulator  is  at  electron  phase  C,  »  cos-1  (-  8  /|a|)—  <(> 

[11].  The  synchroton  oscillation  frequency  for  the  trapped  electrons  is  given  by 

v5  ~  |a|2  -82J'4;  for  the  untapered  where  8  =  0  the  synchrotron  frequency  reduces  to 


[11]. 


The  synchrotron  oscillation  frequency  mixes  with  the  carrier  wave  and  appears  as 


sidebands  around  the  fundamental  at  v0  ±  vs  causing  a  shift  from  the  fundamental 


wavelength  by  A  A/ A  =vs/  2nN.  This  is  often  undesirable  for  FEL  operation,  because 


the  optical  spectrum  broadens.  In  weapons  applications,  the  laser  beam  must  propagate 
through  a  narrow  window  frequency  in  the  atmosphere  where  there  is  small  absorption. 
The  trapped-particle  instability  may  be  removed  by  increasing  the  desynchronism  to 
reduce  the  power  so  that  the  electrons  execute  one  synchrotron  oscillation 
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(V,  =  Va  =  In).  Figure  16  shows  an  example  of  trapped-particle  instability.  Note  that 
the  optical  field  is  strong  a  =  39,  and  a  sideband  appears  on  the  right  of  the  fundamental 
frequency  at  the  optical  spectrum  at  vs=2j. 


-3  2  3  -47  v  47  0  n  1000 


Figure  16.  The  Trapped-Particle  Instability 


G.  LIMIT  CYCLE  BEHAVIOR 

Limit  cycle  behavior  can  occur  in  optical  pulse  structure  and  power  evolution 
when  the  FEL  is  operating  in  strong  optical  fields  with  short  pulses.  In  strong  fields  when 
the  trapped-particle  instability  occurs,  the  oscillations  at  the  synchrotron  frequency  mix 
with  the  optical  carrier  frequency  causing  modulations  of  the  optical  wave  envelope  and 
spectral  sidebands  at  the  synchrotron  frequency.  This  modulation  of  the  optical  pulse 
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shape,  in  combination  with  the  desynchronism  effect,  create  an  optical  pulse  that 
continually  changes  from  one  shape  to  the  other  over  many  passes,  resulting  in  power 
oscillations.  Figure  17  shows  an  example  of  limit  cycle  behavior. 


-8  z  8-47  V  47  0  n  100 

Figure  17.  The  Limit  Cycle  Behavior 

As  it  is  seen  in  the  pulse  evolution  |a(z,n)j,  subpulses  are  continually  being 

formed  at  the  trailing  edge  of  the  optical  pulse  and  march  forward  in  z  over  many  passes 
due  to  the  desynchronism  mechanism.  As  the  subpulses  move  forward  they  enter  a 
region  of  higher  gain,  and  grow,  and  then  continue  their  movement  and  decay  due  to 
resonator  losses.  Limit  cycle  behavior  occurs  at  moderate  values  of  desynchronism  and  is 
usually  undesirable  for  weapons  applications  because  the  optical  spectrum  is  broadened 
by  the  presence  of  the  sidebands  and  because  steady  state  power  is  not  achieved. 
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However,  it  can  easily  be  removed  by  altering  the  desynchronism  as  a  result  of  moving 
one  resonator  mirror  by  microns. 
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III.  THE  TJNAF  FEL  SIMULATIONS 


A.  INTRODUCTION 

Experiments  using  the  TJNAF  FEL  have  been  done  to  explore  the  operation  with 
inversely  tapered  undulators.  In  this  chapter,  numerical  simulations  will  be  presented 
using  the  TJNAF  experimental  parameters.  First,  single  mode  simulations  will  be  used  to 
explore  the  effects  of  negative  taper  on  gain.  Then,  multimode  simulations  will  be  used 
to  analyze  the  operation  of  the  TJNAF  FEL  with  the  negatively  tapered  undulator, 
describing  the  evolution  of  short  optical  pulses  in  the  far  infrared  and  showing  the  effects 
of  taper  and  desynchronism  on  gain  and  power.  In  our  simulations,  the  TJNAF  FEL  will 
be  driven  by  34.5  MeV  and  47.5  MeV  short  electron  pulses.  Tables  2  and  3  show  the 
parameters  for  both  34.5  MeV  and  47.5  MeV  operation. 


37 


Optical  Wavelength  X 

6  pm 

Undulator  Wavelength  Xq 

2.7  cm 

Peak  Undulator  Magnetic  Field  B 

5.5  kG 

Undulator  Periods  N 

41 

Undulator  Length  L 

1.1  m 

Undulator  Parameter  K  (rms) 

0.98 

Taper  Phase  Acceleration  S 

<5  =  0,  -4k,  -6k,  -8ji 

Electron  Energy  ymc2 

34.5  Mev 

Electron  Pulse  Length  az  normalized  to  NX 

1.0 

Initial  Phase  Velocity  vo  for  peak  gain  in  weak  fields 

vo~  2.4  for  3=  0 

vo  =  8.7  for  <5  =  -4jr 

vo  =  1 1 .7  for  <5  =  -6k 

vo~  14.7  for  S  =  -8 k 

Peak  Current  / 

50  A 

Current  Density  j 

10 

Initial  Optical  Field  a0 

0.0 

Cavity  Losses  1/Q 

0.1 

Table  2.  TJNAF  34.5  MeV  Electron  Pulse  Energy  Parameters 
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Optical  Wavelength  A 

3  pm 

Undulator  Wavelength  A0 

2.7  cm 

Peak  Undulator  Magnetic  Field  B 

5.5  kG 

Undulator  Periods  N 

41 

Undulator  Length  L 

1.1  m 

Undulator  Parameter  K  (rms) 

0.98 

Taper  Phase  Acceleration  S 

<5  =  0,  -47t,  -671,  -871 

Electron  Energy  ymc 

47.5  MeV 

Electron  Pulse  Length  az  normalized  to  NX 

1.8 

Initial  Phase  Velocity  v0  for  peak  gain  in  weak  fields 

vo  —  2.4  for  <5  =  0 

vo  =  8.7  for  <5  =  -47t 

vo=  11.7  for  <5  =  -671 

v0=  14.7  for  3  =  -8ti 

Peak  Current  I 

50  A 

Current  Density  j 

7 

Initial  Optical  Field  a0 

0.0 

Cavity  Losses  1/Q 

0.1 

Table  3.  TJNAF  47.5  MeV  Electron  Pulse  Energy  Parameters 
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B.  SINGLE  MODE  SIMULATION  RESULTS 


Single  mode  behavior  was  explored  by  examining  the  effects  of  negative  tapering 
on  gain  for  a  range  of  initial  optical  fields  (0  <  a0  <  40),  and  initial  phase  velocities 

(-30  <  v0  <  30).  The  tapering  phase  acceleration  has  values  8  =  0,-4n-6n  and  -  8;r, 

that  correspond  to  linear  tapering  rate  of  the  dimensionless  undulator  parameter 
A K/K=  0%,  5%,  7.5%  and  10%. 

Figure  18  shows  the  gain  plotted  as  a  function  of  the  initial  optical  field  a0  and 
the  initial  phase  velocities  vo  for  current  j  =  7.  For  3  =  0  and  weak  field  (a0  »  0),  the 
maximum  gain  is  G  =  1 1 3%  and  occurs  at  vo  =  2.4.  For  strong  field  ( a0  =  40),  the 

maximum  gain  reduces  to  G  =  4.7%  and  shifts  more  off  resonance  at  vq=  10.8.  After 
introducing  negative  tapering,  the  maximum  weak  field  gain  reduces  to  G  =  77%  at 
vo=  8.7  for  8  =  -4 k,  G  =  46%  at  vo=  1 1.7  for  3  =  -6k,  and  G  =  19%  at  vo~  14.7  for 
8  =  -8ji.  The  strong  field  gain  increases  to  G  =  1 1%  at  vo  =  1 2.3  for  8  =  -4%,  G  =  9.9%  at 
V0=  15  for  3  =  -6k,  and  G  =  9.5%  at  v0=  24  for  8  =  -8k.  As  the  inverse  taper  increases, 
peak  gain  in  weak  fields  is  shifted  to  the  right  along  the  v0  axis  by  -812  and  decreases, 
whereas  the  strong  field  peak  gain  increases.  As  the  taper  increases,  the  gain  spectrum 
becomes  more  distorted  and  a  second  peak  gain  develops.  For  8  =  -87c,  two  comparable 
peaks  appear,  one  at  vo  —  14.7  and  the  other  at  vo~  23.  Figure  19  shows  the  gain  plotted 
as  a  function  of  the  initial  optical  field  a  and  the  initial  phase  velocities  vo  for  current 
7=10.  The  graphs  are  similar  to  the  j  =  1  curve;  the  only  difference  is  that  the  values  of 
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gain  for/  =  10  are  bigger  by  approximately  52%  in  weak  fields  and  by  approximately 
44%  in  strong  fields. 
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Figure  18.  The  FEL  Gain  Spectrum  GCv^ao)  for/=7 
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Figure  19.  The  FEL  Gain  Spectrum  G(yo,a0)  for/=10 
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C.  MULTIMODE  SIMULATION  RESULTS 

My  first  task  was  to  plot  the  final  steady  state  power  of  the  FEL  as  a  function  of 
desynchronism  d  for  S  =  0,  -4n,  -6n  and  -871.  In  order  to  achieve  this,  we  ran  FEL  pulse 
evolution  simulations  for  various  values  of  desynchronism  d  as  described  in  Chapter  II, 
Section  E.  The  input  parameters  of  the  TJNAF  FEL  operation  used  in  the  simulations  are 
provided  in  Table  2  for  34.5  MeV,  and  in  Table  3  for  47.5  MeV  energy  electron  pulses. 

The  second  task  was  to  plot  the  weak-field  steady-state  gain  as  a  function  of 
desynchronism  d  for  (5  =  0,  -47i,  -671  and  -8tt.  The  same  FEL  pulse  evolution  simulations 
were  run  but  gain  G  was  plotted  as  a  function  of  n  instead  of  plotting  the  power  P.  The 
simulations  were  run  long  enough  to  reach  steady  state  gain  in  weak  optical  fields 
(a  «  K)  with  O  =  1 010  in  order  to  represent  no  losses. 

The  third  task  was  to  plot  the  induced  electron  energy  spread  Ay/y  as  a  function  of 
desynchronism  d  and  taper  phase  acceleration  S.  This  was  accomplished  using  the 
formula  Ay / y  =  Ave  i 4/rN,  where  Ave  is  the  electron  spectrum  full  width  measured  from 

the  electron  distribution  f(v,n)  on  the  FEL  pulse  evolution  simulations,  and  N  =  41  is  the 
number  of  the  undulator  periods. 

Finally,  the  fourth  task  was  to  plot  the  optical  spectrum  width  as  a  function  of 
desynchronism  d  and  taper  phase  acceleration  S,  measured  from  the  power  spectrum 
P(v,n)  on  the  FEL  pulse  evolution  simulations. 
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1.  TJNAF  FEL  Simulation  Results  Using  34.5  MeV  Energy  Electron 
Pulses 

Figure  20  shows  the  graph  of  the  steady-state  power  as  a  function  of 
desynchronism  d  and  taper  phase  acceleration  S.  For  3  =  0,  the  laser  operating  range  is 
from  d  =  0.0035  to  d  =  0.385,  and  reaches  maximum  power  at  d  =  0.004  of  P  =  206.  For 
3  =  -4ji,  the  laser  operating  range  is  from  d  =  0.004  to  d  =  0.29,  and  reaches  the  maximum 
at  d=  0.006  with  P  =  1 89.  For  3  =  -6k,  the  laser  operating  range  is  from  d  =  0.005  to 
d  =  0.19,  and  reaches  the  maximum  at  d=  0.012  with  P  =  148.  Finally  for  3  =  -871,  the 
laser  operating  range  is  from  d=  0.0065  to  d  =  0.053,  and  reaches  the  maximum  at 
d=  0.0125  with?  =  90. 

The  operating  range  decreases  as  the  value  of  3  increases.  It  can  be  seen  that  a 
small  change  in  d  causes  large  difference  in  steady-state  power  for  small  values  of  d.  The 
high  power  area  for  small  d  is  an  area  of  instability.  Increasing  the  value  of  d,  the  FEL 
operation  becomes  more  stable  but  the  power  decreases.  For  the  largest  values  of  d,  the 
power  becomes  even  smaller  due  to  the  reduced  coupling  between  the  optical  and 
electron  pulses. 

The  efficiency  for  the  maximum  values  of  the  steady  state  power  is  r\  =  2.2%  for 
S  =  0,  n  =  1  .9%  for  3  =  -4ji,  r\  =  1 .5%  for  3  =  -6n  and  r\  =  0.94%  for  3  =  -871.  The  highest 
power  is  achieved  from  d  =  0.0035  to  d  =  0.006  and  from  d  =  0.1 8  to  d  —  0.385  for  <5  —  0, 
from  d  =  0.006  to  d  =  0.015  and  d  =  0.04  to  d  =  0.1 8  for  <5  =  -4n;  and  from  d  =  0.015  to 
d  =  0.04  for  3  =  -6n.  Even  though  the  maximum  power  is  obtained  for  3  =  0, 8  =  -4k 
provides  higher  power  at  desynchronism  values  where  the  operation  is  more  stable. 


45 


Figure  20.  Steady-state  Power  vs.  Desynchronism  for  34.5  MeV  Energy 
Electron  Pulses. 
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For  small  values  of  d  where  the  power  is  too  high  and  the  optical  fields  are  too 
strong,  electrons  are  trapped  in  deep  potential  wells  and  perform  synchrotron  oscillations 
causing  the  trapped  particle  instability  is  evident.  For  S  —  0,  this  occurs  between  the 
values  of  d=  0.035  and  d=  0.07,  for  S  =  -4n  between  d  =  0.004  and  d  =  0.02,  for  §  =  -6n 
between  d  =  0.005  and  d=  0.018,  and  for  S  =  -Sn  between  d=  0.009  and  d=  0.022. 
Figures  21  and  22  show  the  pulse  evolution  for  d  =  0.0055  with  the  trapped-particle 
instability  for  <5  =  0  and  S  =  -4n.  By  tapering  the  undulator  the  sideband  gain  is  reduced, 
the  optical  spectrum  width  P(v,n)  becomes  smaller  and  the  steady  state  power  is  also 
reduced. 
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Figure  21.  Trapped  Particle  Instability  Example  for  S  =  0  and  d=  0.0055 
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Figure  22.  Trapped  Particle  Instability  Example  for  d  =  -4n  and  d  =  0.0055 

For  <5  =  0  between  the  values  of  d  =  0.02  and  d  -  0.08,  limit  cycle  behavior  is 
observed.  It  is  also  evident  for  <5  =  -8 k  between  d  =  0.0 1 45  and  d  —  0.015,  and  between 
d  =  0.033  and  d  =  0.C '  Limit  cycle  behavior,  as  described  in  Chapter  II  Section  G,  in 
optical  pulse  structure  and  power  occurs  when  trapped  particles  in  strong  optical  fields 
combine  with  short  optical  pulses.  For  <5  =  0  the  maximum  modulation  is  9.5%  of  the 
average  power  at  d  =  0.04  the  minimum  modulation  is  1 .5%  at  d  =  0.08.  For  S  =  -8ji,  the 
maximum  modulation  is  6.6%  of  the  average  power  at  d=  0.015,  and  the  minimum 
modulation  is  1 .2%  at  d-  0.033.  Black  circles  and  squares  in  Figure  20  indicate  the 
average  value  of  the  power  oscillation.  Figures  23  and  24  show  the  limit  cycle  behavior 


48 


region  for  6  =  0  and  5  =  -8n  while  Figure  25  shows  the  limit-cycle  behavior  pulse 


evolution  for  <5-0  and  d  -  0.04. 


Figure  23.  Limit  Cycle  Behavior  Region  for  S  -  0 


Figure  24.  Limit  Cycle  Behavior  Region  for  S  -  -8jt 
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Figure  25.  Limit  Cycle  Behavior  Pulse  Evolution  for  S  =  0 

At  small  values  of  desynchronism  the  optical  pulse  is  short,  only  one  slippage 
distance  long,  results  in  a  broad  optical  spectrum.  At  large  values  of  desynchronism,  the 
optical  pulse  is  advanced  so  far  each  pass  that  it  decouples  from  the  incoming  electron 
pulse  and  the  steady-state  power  is  reduced.  In  this  region,  the  optical  pulse  is  longer 
than  the  electron  pulse,  sometimes  five  to  six  slippage  distances  long,  and  has  a  long 
exponential  leading  edge  of  the  form  |a(z)|  °c  e~z/4Qd  [10].  The  optical  spectrum  for  large 

desynchronism  and  a  long  pulse  is  narrow.  Figure  26  shows  the  pulse  evolution  for  a 
small  value  of  desynchronism  d  =  0.005  and  S  =  -6%,  and  Figure  27  shows  the  results  for 
large  value  of  desynchronism  d  =  0.16  and  taper  S  —  -6n.  Note  that  the  optical  pulse  in 
Figure  26  is  one  slippage  distance  long  (NX),  and  in  Figure  27  is  six  slippages  distances 
long. 
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Figure  27.  Pulse  Evolution  for  d  =  -6n  and  d  =  0.16 


Figure  28  shows  the  weak-field,  steady-state  gain  as  a  function  of  desynchronism 
d  and  taper  phase  acceleration  S.  The  gain  is  small  for  a  small  d  and  large  d  so  that  many 
passes  are  required  to  reach  steady-state  power.  In  the  middle  of  d  the  gain  is  longer  and 
simulations  require  a  smaller  number  of  passes.  For  <5  =  0,  the  maximum  value  of  the 
weak-field,  steady-state  gain  is  G  =  85%  at  d  -  0.14,  for  <5  =  -4tt  is  G  =  56%  at  d=  0.1 1, 
for  8  =  -6n  is  G  =  32%  at  d=  0.07  and  for  8  =  -8ir  is  G  =  13%  at  d=  0.03. 


Figure  28.  Weak  Field  Steady  State  Gain  vs.  Desyncronism  for  34.5  MeV 
Energy  Electron  Pulses 

For  8  =  -8tt  and  for  values  of  resonator  losses  factor  above  <9=12,  gain 
oscillations  were  observed  between  d  =  0.03  and  d  =  0.09  with  average  modulation 
amplitude  13%  of  the  average  gain.  The  optical  pulse  becomes  spatially  modulated  at  the 
slippage  distance,  but  there  are  no  trapped  electrons  because  the  field  is  not  strong.  This 
is  believed  to  be  a  new  effect  caused  by  the  mode  competition  between  the  two 
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comparable  peaks  in  the  gain  spectrum  for  <5  =  -871.  The  average  values  of  the  gain  in 
limit  cycle  behavior  region  are  indicated  in  Figure  28  with  black  squares.  Figure  29 
shows  the  oscillations  in  gain  for  3  —  -8it  and  d  —  0.045. 


Figure  29.  Weak  Field  Gain  Oscillations  for  6  -  -8n 

Figures  3 1  and  32  show  the  induced  electron  energy  spread  Ay/y,  as  a  function  of 
desynchronism  d  and  tapering  rate  <5.  The  electron  energy  spread  curves  follow  the  trend 

of  the  power  curves.  For  high  power,  there  is  a  large  energy  spread.  For  <5  =  0,  the 

maximum  value  is  8%  at  d  =  0.004,  for  (5  =  -4*  the  maximum  value  is  7.3%  at  d =  0.005, 
for  <5  =  -6n  the  maximum  value  is  6.3%  at  d  =  0.008,  and  for  <5  =  -8tt  the  maximum  value 
is  4.8%  at  d  =  0.0125.  A  goal  of  the  TJNAF  FEL  is  to  keep  Ay/y  less  than  6%  in  order  to 
allow  effective  recirculation  of  the  electron  beam.  The  energy  spread  is  less  than  6%  for 
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desynchronisms  larger  than  d~  0.01  for  all  the  tapering  rates  <5.  The  maximum  power 
attained,  with  energy  spread  less  than  6%,  is  P  =  1 65  for  S  -  - 4 %  at  d  =  0.01. 

Figures  33  and  34  shows  the  optical  spectrum  width  as  a  function  of 
desynchronism  d  and  tapering  rate  <5.  For  <5=0,  the  maximum  width  is  Av  =  35.02  at 
d  =  0.0038,  for  <5  =  -An  the  maximum  width  is  Av  =  24.06  at  d—  0.006,  for  S  =  -6k  the 
maximum  width  is  Av  -  1 5.41  at  d  -  0.008,  and  for  <5  =  -8tt  the  maximum  width  is 
Av  =  16  for  d  =  0.015.  It  can  be  seen  that  tapering  reduces  the  optical  spectrum  width. 

S.  Benson,  J.  Gubeli,  and  G.R.  Neil  conducted  experiments  on  TJNAF  FEL 
oscillator  with  linear  tapering  using  the  same  input  data  as  our  simulations.  The  results  of 
their  experimental  study  are  described  in  Ref  [12].  Figure  30  shows  a  comparison 
between  the  experimental  and  the  simulation’s  FEL  operating  range  in  desynchronism  as 
a  function  of  the  tapering  rate  <5.  The  vertical  axis  plots  the  desynchronism  curve  width  at 
taper  <5  normalized  to  the  desynchronism  curve  width  at  <5  =  0.  As  we  see  the  results 
match  experiment  very  well.  This  is  an  indication  that  the  simulation  results  are  valid. 


Figure  30.  Comparison  Between  Experimental  and  Simulation  Operating 
Ranges 
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Figure  31.  Electrons  Energy  Spread  vs.  Desynchronism  for  34.5  MeV 
Energy  Electron  Pulses 
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Figure  32.  Close-up  of  Electrons  Energy  Spread  vs.  Desynchronism,  from 
d  =  0  to  d  =  0.01  for  34.5  MeV  Energy  Electron  Pulses 
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Figure  34.  Close-up  of  Optical  Spectrum  Width  vs.  Desynchronism,  from 
d  =  0  to  d  =  0.1,  for  34.5  MeV  Energy  Electron  Pulses 
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2  .  TJNAF  FEL  Simulation  Results  Using  47.5  MeV  Energy  Electron 
Pulses 

Figure  35  shows  the  graph  of  the  steady  state  power  as  a  function  of 
desynchronism  d  and  taper  phase  acceleration  <5.  For  <5  =  0,  the  laser  operating  range  is 
from  d  =  0.0038  to  d  =  0.35,  and  reaches  maximum  power  of  P  =  171  at  d =  0.0055.  For  8 
=  -471  the  laser  operating  range  is  from  d  =  0.0044  to  d—  0.27,  and  reaches  maximum 
power  of  P  =  1 80  at  d =  0.007.  For  8  =  -6%,  the  laser  operating  range  is  from  d  =  0.0055 
to  <7=  0.18,  and  reaches  maximum  power  of  P  =  136  at  d  =  0.025.  Finally  for  <5  =  -871,  the 
laser  operating  range  is  from  d=  0.01 1  Xod=  0.055,  and  reaches  maximum  power  of 
P  =  78at  J=0.02. 

The  efficiency  corresponding  to  the  maximum  values  of  steady  state  power  are 
1 .39%  for  8  =  0, 1 .42%  for  8  =  -47t,  1 .05%  for  8  =  -67c,  and  0.64%  for  <5  =  -871.  The 
untapered  undulator  gives  the  best  performance,  for  large  values  of  desynchronism 
d  =  0.16  to  d  =  0.35.  A  negative  taper  rate  of  8  =  -4tc  gives  the  highest  power  and  stability 
in  the  rest  of  the  operating  range  down  to  d  =  0.16,  which  is  a  high  power  area. 

The  trapped-particle  instability  is  evident  again  at  high  power.  For  <5  =  0,  this 
occurs  between  the  values  of  d  =  0.038  and  <7=  0.06,  for  <5  =  -4%  between  d  =  0.0044  and 
d  =  0.02,  and  for  <5  =  -6n  between  d  =  0.0055  and  d  =  0.02.  For  <5  =  -871,  the  maximum 
value  of  the  optical  field  strength  is  a0  =  24.8  and  there  is  no  trapped-particle  instability. 

For  <5  =  0,  between  the  values  of  d=  0.0095  and  d  =  0.065,  limit  cycle  behavior  is 
observed.  The  maximum  modulation  is  7.2%  of  the  average  power  at  d  =  0.01,  and  the 
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minimum  is  0.3%  at  d—  0.06.  Black  circles  in  Figure  35  indicate  the  average  value  of  the 
power  oscillation.  Figure  36  shows  the  region  of  limit  cycle  behavior  for  S  =  0. 


60 


jdMOd  aids  Apeajs 


Figure  35.  Steady-state  Power  vs.  Desynchronism  for  47.5  MeV  Energy 
Electron  Pulses 
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Figure  36.  Limit  Cycle  Behavior  Region  for  S  =  0 

Figure  37  shows  a  comparison  of  relative  width  between  the  experiment  and  the 
simulation  FEL  operating  ranges  as  a  function  of  the  tapering  rate  S.  The  experimental 
data  came  from  Ref  [12]  and  we  see  that  again  the  results  match,  showing  that 
experiments  can  be  simulated  leading  to  reliable  conclusions. 


Figure  37.  Comparison  Between  Experimental  and  Simulation  Operating 
Ranges 
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Figure  38  shows  the  weak-field,  steady-state  gain  as  a  function  of  desynchronism 
d  and  taper  phase  acceleration  S.  For  (5  =  0,  the  maximum  value  of  the  weak-field,  steady- 
state  gain  is  G  =  76%  at  d  =  0.13,  for  S  =  -4k  the  maximum  gain  is  G  =  52%  at  d  =  0.1, 
for  S  =  -6k  the  maximum  gain  is  G  =  3 1  %  at  d =  0.07,  and  for  S  =  -8k  the  maximum  gain 
is  G  =  12%  at  d=  0.036. 


Figure  38.  Weak-field,  Steady-state  Gain  vs.  Desyncronism  for  47.5  MeV 
Energy  Electron  Pulses 

Figures  39  and  40  show  the  fractional  electron  beam  energy  spread,  Ay/y,  as  a 
function  of  desynchronism  d  and  tapering  rate  S.  This  energy  spread  is  induced  by  the 
FEL  interaction.  For  effective  recirculation  of  the  electron  beam  it  is  desirable  to  keep 
Ay/y  less  than  6%.  The  energy  spread  curves  follow  the  same  trend  of  the  power  curves. 
High  optical  power  induces  a  larger  energy  spread.  For  S  =  0  the  maximum  value  is  6.9% 
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at  d  =  0.0049,  for  8  =  -47t  is  6.5 1%  at  d  =  0.0052,  for  8  =  -6n  is  5.5%  at  d  =  0.007,  and  for 
8  =  -871  is  3.97%  at  d-  0.017.  The  energy  spread  is  less  than  6%  for  all  values  of  d for 
8  =  -671  and  8  =  -871.  For  desynchronism  larger  than  d  =  0.007  for  <5  =  0  and  8  =  -4te,  the 
energy  spread  is  also  less  than  6%.  The  maximum  power  attained  with  less  than  6% 
energy  spread  is  P  =  180  for  8  =  -4ti  at  d  =  0.0075. 

Figures  41  and  42  show  the  optical  spectrum  width  Av  as  a  function  of 
desynchronism  d  and  tapering  rate  8.  For  6  =  0,  the  maximum  width  is  Av  =  27  at 
d  =  0.004,  for  8  -  -4ti  the  maximum  width  is  Av  =  18  at  d  =  0.0045,  for  8  =  -6n  the 
maximum  width  is  Av  =  1 3  at  d  =  0.007,  and  for  8  =  -871  the  maximum  width  is  Av  =  8.4 
for  e/  —  0.0 1 1 .  It  can  be  seen  that  tapering  reduces  the  optical  spectrum  width,  and  high 
power  increases  the  width. 
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Figure  40.  Close-up  of  Electrons  Energy  Spread  vs.  Desynchronism,  from 
d  -  0  to  d  =  0.1,  for  47.5  MeV  Energy  Electron  Pulses 
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Figure  41.  Optical  Spectrum  Width  vs.  Desynchronism  for  47.5  MeV  Energy 
Electron  Pulses 
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IV .  LASER  PROPAGATION  IN  ATMOSPHERE  AND  INTERACTION  WITH 

MATTER 

A.  ATMOSPHERIC  COMPONENTS 

The  atmosphere  is  gaseous  and  extends  for  several  hundred  kilometers  above 

Earth.  The  gas  of  the  atmosphere  consists  of  elements  and  compounds.  The  exact 
composition  varies  with  geographic  location  and  altitude  because  both  the  atmospheric 
pressure  and  temperature  change  with  the  vertical  structure  of  the  atmosphere.  The  most 
common  elements  in  the  atmosphere  are  nitrogen  and  oxygen  that  constitute  over  98% 
percent  of  the  atmosphere  by  volume.  Water  vapor,  carbon  dioxide,  nitrous  oxide,  carbon 
monoxide,  and  ozone  are  the  major  radiation  absorbers.  Water  vapor  has  significant 
affect  on  absorption  of  infrared  radiation.  It  is  also  the  most  variable  because  of  the 
evaporation  of  water  from  bodies  of  water  and  condensation  into  clouds  or  dew.  Carbon 
dioxide  does  not  vary  as  much  as  does  water  vapor,  but  tends  to  be  concentrated  around 
large  cities  and  heavy  vegetation  areas.  Carbon  dioxide  is  a  strong  absorber  in  the 
infrared  wavelength  range  3  to  5|xm. 

B.  ATMOSPHERIC  TRANSMISSION 

In  calculating  the  optical  transmission  of  a  laser  radiation  through  the  atmosphere, 

there  are  three  primary  processes  that  affect  the  radiation:  absorption,  scattering  and 
refractive  index  fluctuations  or  turbulence.  The  atmospheric  components  discussed 
earlier  are  related  to  absorption  and  scattering,  while  the  atmosphere’s  variations  of 
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temperature,  pressure  and  density  contribute  to  turbulence.  The  effect  of  these  factors  is  a 
reduction  of  the  power  of  the  laser  irradiation  that  reaches  the  target. 

Absorption  and  scattering  are  usually  grouped  together  under  the  topic  of 
extinction.  Extinction  is  the  attenuation  in  the  amount  of  radiation  passing  through  the 
atmosphere.  Absorption  is  a  process  wherein  a  photon  of  radiation  is  absorbed  by  a 
gaseous  molecule  of  the  atmosphere,  which  translates  to  a  temperature  change.  When 
radiation  is  scattered,  the  direction  of  the  incident  radiation  is  changed  due  to  collision  of 
a  photon  with  an  atmospheric  molecule  or  particle.  This  process  can  be  thought  of  as  the 
atmospheric  molecule  captures  the  incident  radiation  momentarily  and  sends  it  unchanged 
in  all  directions. 

To  understand  how  the  extinction  affects  transmission  of  radiation  through  the 
atmosphere,  think  of  a  single  wavelength  of  incident  radiation  passing  through  an 
infinitesimal  distance  dx.  The  change  in  flux  O  is 

d<$>  =  -jj&dx  (47) 

where  the  incident  radiation  flux  is  represented  by  G>  and  [x  is  the  extinction  coefficient. 
The  negative  sign  indicates  the  reduction  of  power  from  the  absorption  or  scattering 
processes.  In  general,  the  extinction  coefficient  is  made  up  of  two  components: 

fx  =  e  +  £  (48) 

where  s  is  the  absorption  coefficient  and  £  is  the  scattering  coefficient.  Both  the 
absorption  and  scattering  coefficients  depend  on  the  incident  wavelength.  Because  of  the 
interaction  of  the  incident  radiation  with  the  molecules  of  the  atmosphere  the  radiation  at 


70 


the  output  is  changed  to  O  +  JO.  After  integrating  Equation  (1)  it  is  found  that  the 
attenuation  for  a  finite  distance  through  a  homogeneous  medium  is 


<D  =  O0e  *  =  Oor 


(49)  where  T  =  e'^  is  the  transmittance  of  the  atmosphere  over  a  distance  x.  This 
principle  is  known  as  Beer-Lambert  law.  As  it  was  mentioned  above,  the  extinction 
coefficient  is  highly  dependent  on  wavelength,  so  Beer’s  law  is  usually  written  as 

T(A)  =  e~M(X)x  (5°) 


Figure  43  shows  the  transmittance  as  a  function  of  the  wavelength  for  a  horizontal 


path  of  6000  ft  at  sea  level  and  typical  humidity.  As  we  see  the  maximum  transmittance 
occurs  in  the  region  3  to  4.2  pm.  This  explains  why  in  the  next  chapter  the  irradiations 
were  done  at  X  =  3.1pm  where  the  corresponding  transmittance  is  T  =  62%. 


Figure  43.  Atmospheric  Transmittance  for  a  Horizontal  Path  of  6000  ft  at 
Sea  Level  and  Typical  Humidity,  From  Ref.  [5] 
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C.  THERMAL  DIFFUSION  LENGTH 

Suppose  that  a  powerful  laser  that  delivers  sufficient  amount  of  energy  to  the 
target  has  been  built.  Another  parameter  that  has  to  be  considered  is  the  laser’s  spot  size 
at  the  target.  If  this  diameter  of  the  spot  is  smaller  than  the  thermal  diffusion  length  of  the 
material  of  the  target,  the  material  is  able  to  diffuse  the  laser  beam’s  incident  energy  away 
faster  than  the  incoming  energy  can  melt  the  material.  The  thermal  diffusion  length  is 
characteristic  for  each  material  and  determines  its  ability  to  absorb  and  transport  heat  for 
a  given  spot  size. 

Schriempf  gives  a  detailed  derivation  of  the  thermal  diffusion  length  [6].  He 
gives  the  definition  of  the  thermal  diffusion  length,  D,  as  the  distance  required  for  the 
temperature  to  drop  to  1/e  of  the  initial  value.  The  formula  that  gives  the  thermal 
diffusion  length  is 

D  =  (51) 

where 

_  nF^_ AT2 
~  4<D**r 

(52)  is  the  time  required  to  raise  the  material’s  temperature  from  ambient  to  melting  (AT), 
F  is  the  thermal  conductivity,  O0  is  the  initial  radiative  flux  and  k  =  F  /  pFSC  is  a 

constant  where  pFS  is  density  and  C  is  the  specific  heat. 

In  the  next  chapter,  we  are  going  to  describe  the  damage  produced  on  Slip-cast 
Fused  Silica  samples,  so  we  calculate  the  thermal  diffusion  length  for  this  material.  The 
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material  density  is  pFS  =  2200  Kg/m3,  the  specific  heat  is  C  =  920  J/Kg-K,  the  thermal 
conductivity  is  F  =  1 .26  W/m-K,  the  initial  power  density  is  <J>0  =  108  W/m2  and  the 
melting  temperature  of  the  material  is  Tm=  1980  K.  Using  these  values  in  Equation  (5)  it 
is  found  that  the  thermal  diffusion  length  for  Slip-cast  fused  Silica  is  D  —  0.021  mm.  In 
order  to  minimize  thermal  diffusion  and  melt  the  material,  the  laser  spot  size  on  the  target 
must  have  diameter  larger  than  the  thermal  diffusion  length  calculated  above.  This  was 
attained  in  all  the  irradiations  that  were  done  on  the  samples. 

D.  SCALING 

Experiments  indicate  that  the  typical  average  power  density  required  to  kill  an 
incoming  missile  is  about  10  kW/cm2  over  a  spot  area  of  100  cm  .  That  means  that  the 
energy  that  must  reach  the  missile  is  about  P  =  (1  OkW/cm2 )  •  1 00cm2  =  1MW . 
Considering  atmospheric  extinction,  more  energy  is  necessary  at  the  laser  to  obtain  the 
required  energy  at  the  target.  However  the  cost  of  creating  such  a  powerful  laser  is 
extremely  high,  and  in  order  to  come  to  safe  conclusions  about  its  effectiveness,  we  must 
do  several  experiments  using  the  existing  lower  power  FEL. 

The  TJNAF  FEL  is  capable  of  several  hundred  watts  average  power  and 
developments  are  being  made  to  be  upgraded  to  operate  at  10  kW  average  power  in  the 
near  future.  Using  the  existing  power,  we  can  simulate  the  power  density  of  10  kW/cm 
by  focusing  the  beam  to  smaller  spot  sizes,  and  evaluating  the  damage  to  different  sample 
materials.  Thus  in  the  future,  we  are  going  to  be  able  to  develop  scaling  rules  that  will 


73 


allow  us  to  estimate  the  damage  that  a  large  laser  produces  without  spending  enormous 
amount  of  money  building  the  laser  first. 
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V.  FEL  DAMAGE  EXPERIMENTS 


A.  EXPERIMENTAL  PROCEDURE 

Two  experiments  took  place  at  TJNAF  on  slip-cast  fused  silica  samples  on 

August  9, 1999  and  March  14, 2000  by  the  TJNAF  personnel  [1],[2].  In  the  next  few 
paragraphs,  the  experimental  procedures  are  described  the  irradiation  results  are  analyzed. 

In  the  experiment  of  August  9,  1999  two  samples  were  irradiated  through  a 
calcium  fluoride  lens  with  a  measured  back  length  of  137.6  mm.  The  laser  beam 
wavelength  was  2=3.10  pm,  the  pulse  repetition  frequency  (PRF)  was  18.7  MHz  and  the 
power  meter  in  the  optics  control  room  indicated  a  power  of  105  ±  5  W.  The  first  sample 
was  placed  20.7  cm  from  the  back  surface  of  the  lens.  At  this  position,  the  calculated 
waist  radius  of  the  beam  was  0.25  cm,  with  a  corresponding  average  intensity  of  490 
W/cm2.  Three  irradiations  with  no  airflow  were  done,  and  then  the  sample  was  moved  to 
irradiate  fresh  areas,  the  airflow  was  turned  on,  and  three  more  irradiations  were  done. 
The  airflow  was  blowing  across  the  front  face  of  the  sample  at  an  angle  of  0  =  90°  with 
respect  to  the  irradiation  line.  An  Oregon  Scientific  anemometer  was  used,  which 
indicated  a  wind  speed  of  60  mph.  The  irradiation  exposure  time  was  5  seconds. 

After  finishing  the  above  irradiations,  a  new  sample  was  used.  The  sample  was 
moved  in  a  new  position  in  order  to  achieve  a  beam  waist  radius  of  0.087  cm,  which 
yields  to  an  average  intensity  of  10  kW/cm2.  The  same  irradiation  schedule  was 
followed.  The  irradiation  exposure  time  was  again  5  seconds. 


75 


In  the  experiment  of  March  14,  2000,  another  sample  was  irradiated  through  a 
calcium  fluoride  lens  with  a  measured  back  length  of  235.7.  The  laser  beam  wavelength 
was  2  =  3.10  pm,  the  PRF  was  37.4  MHz  and  the  power  meter  in  the  optics  control  room 
indicated  a  power  of  500  ±  10  W.  The  sample  was  placed  217.0  mm  from  the  back 
surface  of  the  lens.  At  this  position,  the  calculated  waist  radius  of  the  beam  was  0.12  cm, 
with  a  corresponding  average  intensity  of  10  kW/cm2.  One  irradiation  with  no  airflow 
was  done,  and  then  the  sample  was  moved  to  fresh  areas,  the  air  was  turned  on,  and  one 
more  irradiation  was  done.  The  air  was  blowing  across  the  front  face  of  the  sample  at  an 
angle  of  0  =  90°  with  respect  to  the  irradiation  line.  An  Oregon  Scientific  anemometer 
was  used,  which  indicated  a  wind  speed  of  83  to  86  mph.  The  irradiation  exposure  time 
was  5  seconds. 

B.  DESCRIPTION  OF  RESULTS 

The  Naval  Research  Laboratory  (NRL)  provided  the  Slip-cast  Fused  Silica 
samples.  The  samples  #1  and  #2  were  new  and  were  used  in  experiments  for  the  first 
time.  The  sample  #3  had  been  used  in  the  past,  and  carried  eight  irradiations  from  a 
previous  experiment. 

1.  Slip-cast  Fused  Silica  Sample  #1 

The  sample  is  2.2  cm  by  7.4  cm  and  has  a  variation  in  thickness  from  0.9  to  1.9 
cm.  It  was  used  for  the  experiments  conducted  on  August  9,  1999.  Figure  44  shows  two 
sets  of  three  irradiations  done  on  the  front  face  of  the  sample.  The  lower  set  of 
irradiations  was  done  while  air  was  blowing  across  the  front  surface  of  the  sample,  and 
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the  upper  set  without  air.  Table  4  shows  the  data  and  the  results  of  these  irradiations. 
Run  number  1,2  and  3  refer  to  the  upper  set  of  irradiations,  and  4,5  and  6  to  the  lower  set 
from  left  to  right. 


Figure  44.  Slip-cast  Fused  Silica  Sample  #1 


Run 

Number 

Wavelength 

(pm) 

PRF 

(MHz) 

Airflow 

(mph) 

Average 

Power 

(Watts) 

Average 

Intensity 

(W/cm2) 

Laser 

Beam 

Diameter 

(mm) 

Damage 

Diameter 

(mm) 

1 

3.10 

18.7 

No 

105 

490 

5 

6.2 

2 

3.10 

18.7 

No 

105 

490 

5 

6.5 

3 

3.10 

18.7 

No 

105 

490 

5 

6.0 

4 

3.10 

18.7 

60 

105 

490  ^ 

5 

5.0 

5 

3.10 

18.7 

60 

105 

490 

5 

5.0 

6 

3.10 

18.7 

60 

105 

490 

5 

5.5 

Table  4.  Irradiation  Data  of  Slip-cast  Fused  Silica  Sample  #1,  After  Ref.  [1] 
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The  exposure  time  of  the  above  irradiations  was  5  seconds  and  as  we  can  see  from 
the  PRF,  the  energy  per  pulse  was  twice  that  of  the  measurements  made  in  March  99. 

The  damage  diameters  were  measured  using  an  optical  microscope.  Due  to  the  low 
average  intensity  (490W/cm2),  we  observed  only  faint  circular  profiles  on  the  sample 
material  after  the  irradiations,  whose  diameters  matched  the  calculated  beam  diameter 
reasonably.  The  damage  was  just  superficial.  The  effect  of  the  airflow  was  to  decrease 
the  diameter  of  the  damage  area.  Figure  45  shows  a  close  up  of  damage  in  Run  2. 
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Figure  45.  Close-up  Damage  to  Slip-cast  Fused  Silica  in  Run  2 


2.  Slip-cast  Fused  Silica  Sample  #2 

The  sample  is  2.2  cm  by  7.4  cm  and  has  a  variation  in  thickness  from  0.9  to  1 .9 
cm.  It  was  used  for  the  experiments  conducted  on  August  9, 1999.  Figure  46  shows  the 
irradiations  done  on  the  front  face  of  the  sample  and  Table  5  the  corresponding  data. 
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Figure  46.  Slip-cast  Fused  Silica  Sample  #2 


Run 

Number 

Wavelength 

(pm) 

PRF 

(MHz) 

Airflow 

(mph) 

Average 

Intensity 

(kW/cm2) 

Laser 

Beam 

Diameter 

(mm) 

Damage 

Diameter 

(mm) 

Penetration 

Rate 

(mm/s) 

1 

3.10 

18.7 

No 

10 

1.76 

3.8 

0.29 

2 

3.10 

18.7 

No 

10 

1.76 

3.5 

0.30 

3 

3.10 

18.7 

No 

10 

1.76 

3.5 

0.32 

4 

3.10 

18.7 

60 

10 

1.76 

3 

0.34 

5 

3.10 

18.7 

60 

10 

1.76 

3 

0.33 

6 

3.10 

18.7 

60 

10 

1.76 

2.5 

0.31 

7 

3.10 

18.7 

No 

10 

1.76 

3.9 

0.29 

Table  5.  Irradiation  Data  of  Slip-cast  Fused  Silica  Sample  #2,  After  Ref.  [1] 


The  exposure  time  for  the  above  irradiations  was  5  seconds  and  the  average  power 
was  105  Watts.  The  average  damage  diameter,  for  the  irradiations  done  in  the  presence 
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of  airflow,  was  2.83  mm  and  for  those  without  airflow  was  3.67  mm.  Thus,  it  can  be  seen 
that  again  the  damage  diameters  matched  the  calculated  beam  diameters  reasonably,  and 
the  effect  of  the  presence  of  air  was  to  decrease  the  diameter  of  the  damaged  area.  No 
bum-through  occurred  during  the  irradiations.  The  damage  produced  on  the  sample  had 
the  shape  of  small  circular  crater.  The  areas  around  the  craters  were  clean  of  debris.  Part 
of  the  melted  material  was  evaporated  during  the  irradiation,  and  the  rest  of  it  remained 
inside  the  crater. 

The  average  penetration  rate,  for  the  irradiations  done  in  the  presence  of  airflow, 
was  0.326  mm/s  and  for  those  without  airflow  was  0.3  mm/s.  The  presence  of  airflow 
slightly  increased  the  penetration  rate. 

Figure  47  shows  the  effects  of  the  exposure  time  on  penetration  depth  rate,  of  an 
older  experiment  conducted  on  a  sample  of  Slip-cast  Fused  Silica. 
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Figure  47.  Exposure  Time  vs.  Penetration  Rate  for  Fused  Silica,  From 
Ref.  [3] 
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The  average  power  of  that  experiment  was  100  to  103W  ±  5W,  the  wavelength 
X  —  4.825  pm,  the  average  intensity  10  kW/cm2  and  the  PRF  37.4  MHz.  As  it  can  be  seen 
from  Figure  4,  the  penetration  depth  rate  that  corresponds  to  exposure  time  5  seconds  is 
0.26  mm/s.  Therefore  when  the  PRF  is  18.7  MHz  and  the  wavelength  is  X  =  3.10  pm,  the 
penetration  rate  is  bigger  either  with  the  presence  of  airflow  or  without  it.  It  is  true  that 
when  the  PRF  is  lower  there  is  higher  fluence  per  pulse  and  thus  intensity  per  micropulse, 
but  we  can’t  tell  for  sure  if  this  improvement  in  the  penetration  rate  is  due  to  the  different 
wavelength  or  the  PRF. 

In  the  future,  it  might  be  interesting  to  conduct  experiments  changing  just  one 
parameter  in  each  experiment  so  that  we  can  evaluate  better  the  effect  of  each  of  them.  It 
would  also  help  in  scaling,  if  in  each  experiment  we  plot  curves  of  exposure  time  vs. 
penetration  depth  rate.  Figure  48  shows  a  close-up  of  damage  in  Run  4. 


Figure  48.  Close-up  of  Damage  of  Slip-cast  Fused  Silica  in  Run  4 
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3.  Slip-cast  Fused  Silica  Sample  #3 

The  sample  is  6.9  cm  by  7.4  cm  and  has  a  variation  in  thickness  from  0.9  to  1 .9 
cm.  It  was  used  for  the  experiments  conducted  on  March  14,  2000.  Figure  49  shows  the 
irradiations  done  on  the  front  face  of  the  sample  and  Table  6  the  corresponding  data. 


Figure  49.  Slip-cast  Fused  Silica  Sample  #3 


Run 

Number 

Wavelength 

(pm) 

P.R.F. 

(MHz) 

Airflow 

(mph) 

Average 

Intensity 

(kW/cm2) 

Laser 

Beam 

Diameter 

(mm) 

Damage 

Diameter 

(mm) 

Penetration 

Rate 

(mm/s) 

1 

3.10 

37.425 

83-86 

10 

2.4 

1.32 

2 

3.10 

37.425 

No 

10 

2.4 

5.6 

7.5 

Table  6.  Irradiation  Data  of  Slip-cast  Fused  Silica  Sample  #3,  After  Ref.  [2] 
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The  average  power  was  500  Watts  and  bum-through  of  the  material  occurred 
during  the  above  irradiations.  The  bum-through  time  for  Run  1  was  6.8  seconds  and  for 
Run  2  was  1 .2  seconds.  In  Run  1 ,  the  damage  diameter  on  the  front  face  of  the  material 
sample  was  4.4  mm,  and  on  the  back  face  was  3.1  mm.  In  Run  2,  the  damage  diameter 
on  the  front  face  was  5.6  mm  and  on  the  back  face  1.8  mm.  It  can  be  seen  that  the  effect 
of  airflow  was  to  reduce  the  front  face  damage  diameter  and  to  increase  the  bum-through 
time. 

It  can  also  be  seen  that  the  back  face  damage  diameters  of  the  material  sample  are 
smaller  than  the  front  face.  There  are  three  reasons  that  explain  why  this  happens.  The 
first  reason  is  that  the  beam  profile  follows  the  Gaussian  distribution,  with  the  highest 
intensity  in  the  center  of  the  beam  and  intensity  down  by  1/e  at  the  beam  radius.  The 
second  reason  is  the  position  of  the  sample.  The  sample  is  9  mm  thick  and  was  located 
during  the  irradiations  217.0  mm  from  the  back  surface  of  the  lens  while  the  back  focal 
length  of  the  lens  was  235.7  mm.  It  was  1 8.7  mm  from  the  focus,  so  the  beam  size  is 
decreasing  as  it  proceeds  through  the  material.  The  third  reason  is  that  the  front  surface  is 
exposed  to  irradiation  for  longer  time  than  the  back. 

Since  bum  through  occurred,  we  have  the  opportunity  to  have  a  better  view  of  the 
damage  and  extract  some  extra  data  about  it.  The  volume  of  the  total  damaged  region  in 
Run  1  (labeled  IT)  is  estimated  by 

9mm 

VlT  =  \n-R\T{z)-dz,  (53) 

0 

where  z  is  the  material  thickness  and  the  radius  changes  approximately  linearly  as 
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Rn  (z)  =  2.2mm  -  0.072  •  z[mm]. 


(54) 


The  volume  of  the  total  damaged  region  in  Run  2  (labeled  27)  is  estimated  by 

9  mm 

V2T  =  \n-R22T(z)-dz,  (55) 

o 

where  the  radius  changes  approximately  linearly  as 

R1T(z)  =  2.8mm  -0.211  -2[mm].  (56) 

When  we  say  entire  damaged  region,  we  mean  both  the  hole  and  the  melted  and 
rehardened  portion.  After  doing  the  above  calculations  the  volume  of  the  entire  damage 
region  in  Run  1  is  Vn=  100  mm  and  in  Run  2  is  V2t~  105  mm  .  Knowing  that  the 
density  of  the  fused  silica  is  pFS  =  2.2  gm/cm3,  the  mass  of  the  entire  damaged  region  in 

Run  1  is  miT~  0.221  gm  and  in  Run  2  is  m2T=  0.232  gm. 

The  volume  of  the  hole  in  Run  1  (labeled  1H)  is  estimated  by 

9  mm 

Vw  =  \zc-R2H{z)-dz,  (57) 

0 

where  the  radius  changes  approximately  linearly  as 

7?1//(z)  =  1.6mm-0.1167-z[mm].  (58) 

The  volume  of  the  hole  in  Run  2  (labeled  2H)  is  estimated  by 

9  mm 

V2H  =  j 'ir-R22H{z)-dz,  (59) 

o 

where  the  radius  changes  approximately  linearly  as 

R2h  (z)  =  0.75mm  -  0.0444  •  z[mm],  (60) 

From  the  above  calculations  the  volume  of  the  hole  in  Run  1  is  Vjh  =  35  mm  and 
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in  Run  2  is  V2H  =  8.9  mm3.  The  density  of  the  fused  silica  is  pFS  =  2.2  gm/cm3,  so  the 
mass  of  the  material  removed  creating  a  hole  in  Run  1  is  miH=  0.078  gm  and  in  Run  2  is 
m2H~  0.020  gm. 

What  it  is  seen  from  these  calculations  is  that  the  damaged  regions,  either  with  or 
without  the  presence  of  airflow,  have  approximately  the  same  volume  (VJT  =100  mm 
and  V2t  =  105  mm3).  The  basic  effect  of  the  airflow  is  that  it  increases  the  volume  of  the 
hole  of  the  damaged  region  (Vih=  35  mm  and  V2h=  8.9  mm  ). 


Table  7  shows  the  bum-through  irradiation  data  of  an  older  experiment  conducted 
on  the  same  sample  of  fused  silica  without  the  presence  of  air  and  analyzed  in  Ref  [3]. 


Run 

Number 

Wavelength 

(pm) 

PRF 

(MHz) 

Average 

Power 

(Watts) 

Average 

Intensity 

(kW/cm2) 

Volume 

of  entire 

damaged 

region 

(mm  ) 

Volume 

of  the 

hole 

(mm3) 

Penetration 

Rate 

(mm/s) 

2 

4.825 

37.4 

100 

10 

92 

5.6 

0.081 

Table  7.  Irradiation  Data  of  Slip-cast  Fused  Silica  Sample  #3,  After  Ref.  [3] 

Comparing  our  irradiation  results  with  those  on  Table  7,  we  see  that  shifting  the 
wavelength  from  A  =  4.825  pm  to  A  =  3.1  pm  and  increasing  the  average  power  from  100 
Watts  to  500  Watts,  the  penetration  rate  increases  from  0.081  mm/s  to  1.32  mm/s  with  the 
presence  of  airflow  and  to  7.5  mm/s  without  it.  The  volume  of  the  entire  damaged  region 
increases  from  92  mm3  to  100  mm3  with  the  presence  of  airflow  and  to  105  mm3  without 
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it.  Finally  the  volume  of  the  hole  increases  from  5.6  mm3  to  35  mm3  with  the  presence  of 
airflow  and  8.9  mm3  without  it. 

Figure  50  shows  a  close-up  of  damage  in  Run  1  and  Figure  51  in  Run  2.  Figure 
52  shows  the  view  of  those  irradiations  on  the  back  face  of  the  sample. 


1mm 


Figure  50.  Close-up  of  Damage  of  Slip-cast  Fused  Silica  in  Run  1 


Figure  51.  Close-up  of  Damage  of  Slip-cast  Fused  Silica  in  Run  2 
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Figure  52.  Damage  Back  View  of  Slip-cast  Fused  Silica  in  Run  1  and  2 
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VI.  CONCLUSIONS 


Anti-ship  missiles  (ASM)  constitute  a  major  threat  to  warships.  Advanced 
technology  is  used  in  the  missile  production,  making  them  variously  faster,  stealthier, 
lower  flying,  less  susceptible  to  countermeasures  and  more  agile  in  the  terminal  phase. 

As  a  result,  a  highly  sophisticated  close  in  weapon  systems  (CIWS)  is  needed.  The 
Phalanx  gun  is  the  most  widely  deployed  CIWS.  Through  a  series  of  Phalanx  simulations 
we  showed  that  it  is  inadequate;  under  ideal  conditions,  the  distance  where  a  hard-kill  is 
achieved  is  too  close  to  the  ship.  Even  though  the  ship  may  not  be  hit  by  an  intact 
missile,  it  cannot  avoid  damage  by  the  missile  debris. 

High-energy  lasers  (HEL)  might  be  an  answer  to  the  highly  sophisticated  ASM, 
and  the  FEL  is  a  candidate  for  CIWS.  The  advantages  that  it  provides  are  almost 
instantaneous  reaction  at  the  speed  of  light,  a  large  missile  destruction  range,  rapid  re¬ 
engagement,  a  large  and  renewable  magazine,  line  of  sight  accuracy,  precision  aim 
pointing,  high  single  shot  cost-effectiveness,  and  tunability  to  specific  optimal 
wavelengths.  The  tunability  is  the  main  advantage  over  other  HELs,  because  a 
wavelength  can  be  selected  that  would  give  good  atmospheric  transmission  or  is 
appropriate  to  the  specific  target  absorption  characteristics. 

The  Navy’s  Directed  Energy  Office  is  currently  funding  TJNAF  to  study  the  FEL 
as  a  possible  ship-defense  weapon.  The  TJNAF  FEL  is  capable  of  several  hundred  watts 
average  power  and  new  developments  are  being  made  to  upgrade  to  operate  10  kW 
average  power  in  the  near  future.  This  thesis  also  studied  the  laser  damage  produced  by 
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TJNAF  FEL  on  Slip-cast  Fused  Silica  samples;  a  typical  material  for  missile  construction. 
Using  the  existing  FEL,  we  estimated  the  damage  from  a  MW-class  weapon  by  focusing  a 
lower  power  beam  to  a  smaller  spot  size.  Comparing  our  experimental  results  with 
previous  experiments  on  the  same  materials,  we  found  that  a  change  in  wavelength  from 
X  =  4.825  pm  to  X  -  3.10  pm  and  changing  the  pulse  repetition  frequency  (PRF)  from  37.4 
MHz  to  18.7  MHz  gave  an  improvement  in  penetration  rate  of  15.4%  without  airflow  and 
25.4%  with  airflow.  In  both  experiments,  the  average  power  was  lOOWatts  focused  to  an 
intensity  of  the  irradiation  was  10  kW/cm2.  Changing  the  wavelength  from  X  =  4.825  pm 
to  X  =  3.10  pm,  increasing  the  power  from  100  Watts  to  500  Watts,  and  keeping  the  PRF 
constant  at  37.4  MHz  with  10  kW/cm2  irradiation  intensity,  the  improvement  in  the 
penetration  rate  was  1530%  with  the  presence  of  airflow  and  9200%  without  airflow. 

The  volume  of  total  damage  increases  by  8.7%  with  the  presence  of  airflow  and  by  14.1% 
without  it,  and  the  volume  of  the  damage  hole  is  increased  by  525%  with  the  presence  of 
airflow  and  by  59%  without  it. 

In  the  future,  it  might  be  interesting  to  conduct  experiments  changing  just  one 
parameter  in  each  experiment  so  that  we  can  better  evaluate  the  effect  of  each  of  them. 

The  eventual  goal  of  the  damage  experiments  is  to  develop  scaling  rules  that  will  reliably 
predict  the  damage  from  a  larger  laser  without  spending  enormous  amount  of  money 
building  that  laser  first. 

Experiments  using  the  TJNAF  FEL  have  explored  the  operation  with  inversely 
tapered  undulators.  In  this  thesis,  we  also  described  the  single  mode  and  multimode 
simulation  results,  using  the  TJNAF  FEL  experimental  parameters.  Steady-state  power. 
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weak-field  steady-state  gain,  induced  electron  energy  spread,  and  optical  spectrum  width 
behaviors  were  explored  as  a  function  of  desynchronism  d  and  tapering  rate  8,  using  FEL 
pulse  evolution  simulations  including  short  pulse  effects. 

For  FEL  operation  driven  by  the  34.5  MeV  energy  electron  pulses,  it  was  found 
that  even  though  the  maximum  power  is  obtained  for  8  =  0,  a  negative  taper  of  8  =  -4n 
provides  higher  power  that  is  more  stable.  The  maximum  power  attainable  while  keeping 
induced  energy  spread  less  than  6%  is  P  =  1 65  for  negative  taper  of  <5  =  -4n  at  d  =  0.01 . 

For  the  FEL  operation  driven  by  the  47.5  MeV  energy  electron  pulses,  it  was 
found  that  negative  tapering  with  tapering  rate  8  =  -4iz  gives  the  highest  power  and 
stability  up  to  d  =  0. 1 6.  The  maximum  power  attainable  while  keeping  the  induced 
energy  spread  less  than  6%  is  P  =  177  for  8  =  -4n  at  d=  0.0075.  Keeping  Ay/y  less  than 
6%  is  the  requirement  for  effective  recirculation  of  the  electron  beam. 
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